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ABSTRACT 

As part of an extensive study of the physical properties of active galactic nuclei 
(AGN) we report high spatial resolution near-IR integral-field spectroscopy of the 
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narrow-line region (NLR) and coronal-line region (CLR) of seven Seyfert galax- 
ies. These measurements elucidate for the first time the two-dimensional spatial 
distribution and kinematics of the recombination line Br7 and high-ionization 
lines [Si Vl], [Al ix] and [Ca vm] on scales < 300 pc from the AGN. The observa- 
tions reveal kinematic signatures of rotation and outflow in the NLR and CLR. 
The spatially resolved kinematics can be modeled as a combination of an outflow 
bicone and a rotating disk coincident with the molecular gas. High-excitation 
emission is seen in both components, suggesting it is leaking out of a clumpy 
torus. While NGC 1068 (Seyfert 2) is viewed nearly edge-on, intermediate-type 
Seyferts are viewed at intermediate angles, consistent with unified schemes. A 
correlation between the outflow velocity and the molecular gas mass in r < 30 
pc indicates that the accumulation of gas around the AGN increases the colli- 
mation and velocity of the outflow. The outflow rate is 2-3 orders of magnitude 
greater than the accretion rate, implying that the outflow is mass loaded by the 
surrounding interstellar medium (ISM). In half of the observed AGN the kinetic 
power of the outflow is of the order of the power required by two-stage feedback 
models to be thermally coupled to the ISM and match the Mbh — °"* relation. 
In these objects the radio jet is clearly interacting with the ISM, indicative of a 
link between jet power and outflow power. 

Subject headings: galaxies: active - galaxies: Seyfert - galaxies: nuclei - galaxies: 
kinematics and dynamics- line: profiles - infrared: galaxies 



1. Introduction 

The narrow-line region (NLR) in active galactic nuclei (AGN) plays a crucial role in 
attempts to understand the AGN structure and evolution, as well as the interaction between 
the nuclear engine and the circumnuclear interstellar medium (ISM) of the host galaxies. 
Among the fundamental AGN components, the NLR is the largest observable structure in 
the UV, optical and near-IR that is directly affected by the ionizing radiation and dynamical 
forces from the active nucleus. The size and properties of the NLR has been scrutinized 
intensely over the past few decades. Ground-based and Hubble Space Telescope (HST) 
[O ill] A5007A (hereafter [O ill]) and Ha + [N n] narrow-band imaging of Seyfert galaxies 
(Pogge 1988a, 1988b; Haniff et al. 1988; Mulchaey et al. 1996; Evans et al. 1991, 1993; 
Schmitt et al. 1996, 2003) revealed extended emission in several of these objects. The overall 
morphology of the NLR is varied, but there is a trend in the sense that the NLRs in Seyfert 2 
galaxies are more elongated than those in Seyfert 1 galaxies. The fact that several Seyfert 2s 
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presented conically shaped NLRs (e.g., NGC 1068 and NGC 4388) indicates that the source 
of radiation ionizing the gas is collimated, a result that gave strong support to the torus 
version of the Unified Model of AGN (Antonucci 1993; Urry & Padovani 1995). However, 
when comparing ground-based [O ill] images of Seyfert galaxies from Mulchaey et al. (1996) 
with the HST snapshot survey of Schmitt et al. (2003), the latter reveal, on average, four 
times smaller NLR sizes, probably due to low sensitivity (quick exposures). Furthermore, 
one needs to be cautious when talking about "NLR size" as [O ill], or other low- ionization 
lines, can be produced by other physical processes not related to the AGN such as star 
formation or shock- ionized gas. 

The physical conditions in the NLR have also been studied extensively, and it is now 
widely accepted that photoionization by the central continuum source plays a dominant role 
in exciting the gas (e.g. Kraemer & Crenshaw 2000a; Kraemer et al. 2000; Bennert et al. 
2006) although shocks may play an important role in localized regions (Kraemer & Crenshaw 
2000b; Dopita et al. 2002; Contini et al. 2009). 

Ground-based long-slit spectroscopic studies attempted to determine the kinematics of 
the NLR, but a general consensus on the velocity flow pattern was not reached; cases were 
made for infall, rotation, outflow, etc. (e.g. Veilleux 1991; Moore & Cohen 1996; Cooke et 
al. 2000; Fraquelli et al. 2000; Ruiz et al. 2001; Fraquelli et al. 2003). Spectroscopic studies 
with the Faint Object Camera (FOC) on board of HST also fail in this regard. Several 
models were advocated on the basis of these data, including rotation, radial acceleration by 
radio jets, and lateral expansion of gas around the jets (Winge et al. 1997, 1999; Axon et al. 
1998; Capetti et al. 1999). Observations with the Space Telescope Imaging Spectrograph 
(STIS) on HST of a number of key objects have suggested that radial outflow is the dominant 
component of motion in the NLR. The first reliable evidence for outflow in the NLR came 
from HST /STIS long slit observations of NGC 1068 (Crenshaw & Kraemer 2000) and NGC 
4151 (Hutchings et al. 1998; Crenshaw et al. 2000). Since then, detailed constraints on 
kinematic models of the NLR in Seyfert galaxies became possible (see, e.g., Crenshaw et al. 
2010 and references therein). Although the comparison of the radial velocity curves with 
biconical models provide a good fit to the overall flow pattern, there exist significant local 
variations that remain unexplained. Furthermore, the dynamics, origin and energetics of the 
outflows, as well as the role of the radio jet and the presence of other types of motion in the 
NLR kinematics, are not fully understood. Two-dimensional velocity fields with high spatial 
resolution are needed to get a more complete picture of the kinematics in the NLR. Recent 
studies based on integral-field spectroscopy (a few of them done with Adaptive Optics in 
the near-IR) have started to find evidence for radial outflows in the NLR of Seyfert galaxies 
(Garcia-Lorenzo et al. 2001; Barbosa et al. 2009; Storchi-Bergmann et al. 2010; Riffel et al. 
2010). 
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The high-ionization lines, or coronal lines, observed in the spectra of many AGNs, are 
forbidden fine structure transitions in the ground level of highly ionized atoms which have 
threshold energies above 100 eV. In contrast to [O ill] or Ha, which are commonly used 
to study the NLR, the coronal lines are free of potential contributions from star formation. 
Therefore, these lines provide direct insights into the structures and dynamical forces as- 
sociated with the AGN. Ground-based spectroscopic studies revealed that these lines are 
present with comparable strength in AGN spectra regardless of their class (Penston et al. 
1984; Erkens et al. 1997; Prieto & Viegas 2000; Rodriguez-Ardila et al. 2002; Reunanen 
et al. 2003), tend to be broader than NLR lines (400 <FWHM< 1000 km s" 1 ), and their 
centroid velocity is often blueshifted with respect to the systemic velocity of the galaxy (e.g. 
Rodriguez-Ardila et al. 2006). These facts have been interpreted as evidence for a location 
of the coronal-line region (CLR) closer to the AGN than the NLR, but outside the broad 
line region (BLR), and probably associated with outflows (e.g. Miiller- Sanchez et al. 2006, 
hereafter MS06). Because of their high ionization potential (IP> 100 eV), the lines can 
be produced by either a hard UV to soft X-ray continuum, fast shocks, or a combination 
of both. The success of photoionisation models to reproduce, within a factor of few, the 
measured line ratios points towards photoionization as the main excitation mechanism of 



the C LR (IKraemer et all 1200(1 iGroves et all 12004! : Mullanev et al. I I2009J : Mazzalav et al. 



2010 ). although an additional contribution from fast shocks must be invoked to fully explain 
the line ratios in several cases (MS06; Rodriguez-Ardila et al. 2006; Contini et al. 2009). 
However, the morphology, size and kinematics of the CLR remain uncertain because in most 
observational studies the coronal lines are not spatially resolved. 

Detailed studies of the physical conditions across the CLR were not possible bef o re the 
advent of the HST and ground-based Adaptive Optics (AO). Recently, iPrieto at al.l (120051 ) 
determined, for the first time, the size and morphology of the CLR in four Seyfert galaxies. 
These authors found that the CLR as traced by the [Si vn] A2.48/im line extends out to 
r = 30 — 200 pc, and is aligned preferentially with the direction of the lower ionization 
cones seen in [O ill] or Ha. The morphology and size of the CLR indicate that the nuclear 
radiation is either intrinsically collimated or that some additional excitation mechanism (such 
as fast shocks) contributes to the formation of coronal gas at those o bserved distances from 
the ionizing source. With the long-slit capability of HST (STIS), iMazzalay et al.l (120 10f ) 
attempted to study the CLR kinematics in a sample of ten Seyfert galaxies. They found 
that the radial velocity curves of the coronal lines follow the same pattern as those of [O ill] 
emission and concluded that neither the strength nor the kinematics of the CLR scale in any 
obvious and strong way with the radio jet. 

There is great current interest in NLR and CLR kinematics, to resolve the much debated 
question of feedback from AGN via outflows of ionized gas, now recognized as a potentially 



-5- 



crucial input to the Mbh — o"* relation (e.g. Gebhardt et al. 2000): AGN outflows provide 
energy into the ISM of the host galaxy, controlling the star formation and accretion processes, 
thus regulating the black hole and bulge growth. Recent results from Krongold et al. (2007) 
and Storchi-Bergmann et al. (2010) indicate that the total kinetic energy released by the 
AGN outflow is comparable to the energy required to disrupt the hot phase of the ISM, 
which could then disrupt star formation. These important topics had not previously been 
widely studied because high signal-to-noise (S/N) spectra over a full two-dimensional (2D) 
field with high spatial resolution are required, and these are still scarce. Particularly, as the 
coronal lines are exclusively driven by the AGN power, the CLR is an ideal laboratory for 
investigating the physical processes within < 300 pc from the central black hole. 

In an effort to draw general conclusions about the kinematics of the NLR and CLR in 
active galaxies, we have carried out a program of high spatial resolution, near-IR integral 
field spectroscopy of Seyfert galaxies. The measurements discussed here were obtained with 
the use of AO and the instruments SINFONI on the European Southern Observatory (ESO) 
8.2 m Very Large Telescope (VLT) and OSIRIS on the 10 m Keck II Telescope. In this paper, 
which is the first of a series devoted to investigate the role of NLR and CLR kinematics in 
the co-evolution of supermassive black holes (SMBHs) and their host galaxies, we discuss the 
NLR and CLR distribution and kinematics of seven Seyfert galaxies using the Br7 A2.16/im 
recombination line as tracer of the NLR and the near-IR coronal lines: [Si vi] A1.96/iin (IP 
= 167 eV), [Al ix] A2.04/!m (IP = 285 eV) and [Ca vm] A2.32 ; um (IP = 127 eV). The existing 
studies of the kinematics of the NLR are mostly based on long-slit spectroscopy obtained in 
the optical with HST. Our goal is to extend these studies to the near-IR, a region which 
is less affected by dust, and where AO permits diffraction-limited integral field spectroscopy 
of various emission lines that are a vital addition to the lines observed at optical and UV 
frequencies. Furthermore, this work serves as the basis for further studies of the NLR with 
future near-IR telescopes and interferometers, such as, e.g. LBT, ELT and JWST. We 
aim to determine for the first time the 2D kinematics of the CLR, examine previous claims 
of radial outflow along the NLR, as well as investigate possible interactions between the 
NLR/CLR and other components of the AGN or the host galaxy in the central < 300pc. 
When a significant outflowing structure is present, the 2D kinematics are characterized with 
appropriate kinematic models and the mass outflow rates are calculated. 

The paper is structured as follows. In Section 2 we describe the sample characteristics, 
observations and the methods used for the data reduction and analysis. The general prop- 
erties of the NLR and CLR are described in Section 3. In Section 4 we analyze in detail the 
observed kinematics and perform kinematic modeling. The implications of the results are 
discussed in Section 5, including the origin of the ionized gas motions and possible interac- 
tions with other components of the unified model of AGN and the host galaxy. The overall 
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conclusions of the study are outlined in Section 6, and in the Appendix a more detailed 
description of the individual objects is presented. Throughout this paper, H = 75 km s" 1 
Mpc^ 1 is assumed. 



2. Sample, Observations and Data Processing 

2.1. Galaxy Sample 

The seven galaxies studied in this paper were observed with AO and the instruments 
VLT/SINFONI and Keck/OSIRIS. Source selection of the SINFONI and OSIRIS subsamples 
was driven principally by technical considerations, being the primary criteria for selecting 
AGN that (1) the nucleus should be bright and compact enough for good AO correction, 
used either as natural guide star (NGS) or as a tit-tilt guide for the laser guide star (LGS) 
AO system, (2) the galaxy should be close enough (D < 70 Mpc) that small spatial scales 
can be resolved at the near-IR diffraction limit of an 8-m class telescope, and (3) the galaxies 
should have been observed at different wavelengths so that a multiwavelength approach can 
be used when analyzing the integral field data. The only exception is NGC 6814 in the 
SINFONI subsample which was observed without the AO module. 

The resulting sample of seven galaxies is listed in Table [TJ Both the OSIRIS and SIN- 
FONI subsamples contain NGC 6814 and NGC 7469. The measured and derived properties 
from both datasets are consistent, and a comparison of the data from each instrument is pre- 
sented in the Appendix. Although the OSIRIS data have a better spatial resolution than the 
SINFONI data of these two galaxies, we decided to perform the kinematic analysis using the 
SINFONI datacubes, since the SINFONI FOV covers better the extension of the NLR and 
CLR (see Section [3]), making the fitting procedure more accurate (Section Hj). As can be seen 
in Table [Tj there exists bias towards intermediate type 1 objects. This bias is unintenional, 
as there is no preference for a particular type of Seyfert nucleus for the study presented 
here. Narrow Br7 and coronal lines have been detected in the optical and infrared spectra 
of type 1 and type 2 Seyfert galaxies in approximately equal fractions (Rodriguez- Ardila et 
al. 2011). As can be seen in Table [H we detect broad Br7 (FWHM > 1000 km s" 1 ) in all 
galaxies classified as type 1.5 and 1.9. In general, we only consider an emission-line detected 
and "useful" for kinematic analysis when its power is a factor of 3 above the rms scatter of 
the spectrum (S/N > 3). Narrow Br7 and [Si vi] fulfill this requirement in all galaxies of 
the sample, and therefore will be used as tracers of the NLR and CLR, respectively. 
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2.2. Observations and Data Reduction 

SINFONI and OSIRIS are AO-assisted, cryogenic near-IR integral-field spectrographs, 
commissioned at the VLT and the W. M. Keck Observatory, respectively. They deliver 
spectra simultaneously over a contiguous two-dimensional field of view in the wavelength 
range from 1.1 — 2.45/xm at a resolving power of 1500 — 4000. Details about the SINFONI 
instrument can be found in Eisenhauer et al. (2003) and Bonnet et al. (2004) , and a more 



compl ete discussion on the integral field spectrometer OSIRIS can be found in lLarkin et al. 



(|2006l). 



// 



The SINFONI data presented here were taken in nine separate runs between July 2004 
and June 2009. The spectral range was similar for all the runs (1.95-2.45 //m) with typical 
spectral resolution R~4000 (FWHM 70 km s" 1 ). The pixel scale was either 0.0125" x 0.025 
or 0.05" x 0.1", resulting in a field of view (FOV) of 0.8" x 0.8" or 3.2" x 3.2", respectively. 
The only exception is NGC 6814 which was observed without the AO module and the largest 
pixel scale of 0.125" x 0.25" (FOV=8" x 8"). For each gala xy, the nucleus was used as NGS for 



a near-diffraction-limited correction by the AO module ( iBonnet et al.l 120031 ). The exposure 
times for Circinus, NGC 3783, NGC 1068 and NGC 6814 were 30, 40, 110, and 140 minutes, 
respectively, and for the remaining galaxies 60 minutes of on-source data were obtained. 
Data reduction wa s accomplished using the SINFONI custom reduction package SPRED 



( lAbuter et al.ll2005l ). This performs all the usual steps needed to reduce near-IR spectra, but 
with the additional routines for reconstructing the data cube. The reduction steps include 
trimming, bias and background subtraction, flat fielding, cosmic rays cleaning, correction for 
nonlinearity, alignment and interpolation of the data, extraction of the spectra, wavelength 
calibration, sky subtraction and coaddition of different exposures. Telluric correction and 
flux calibration were performed using a standard star (A- or B-type) observed next in time 
to the science frames. In addition, flux calibration was cross-checked with VLT/NACO data 
in 1" — 3" apertures. In some objects the flux calibration was further cross-checked in larger 
3" — 5" apertures using 2MASS data. Agreement between the different data sources was 
consistent to 20%. 

The OSIRIS galaxies were observed between April 2006 and January 2008 using both 
the NGS and LGS AO systems (van Dam et al. 2004; van Dam et al. 2006; Wizinowich et 
al. 2006). The data have 0.035" pixel" 1 , a FOV of 0.56" x 2.24" and a spectral resolution of 
R~ 3000 in the wavelength range 1.965 — 2.381/mi. Details on the observations of individual 
objects can be found in Hicks et al. (2009, hereafter H09). The OSIRIS data was reduced 
using the OSIRIS final data reduction pipeline (DRP), which performs all the usual steps 
needed to reduce near-IR spectra, but with the additional routines for reconstructing the 
data cube (see above). 



The spatial resolution of each data cube was estimated from both the broad Br7 emission 
and the non-stellar continuum, the two always unresolved in 8-m class telescopes. The 
only exception was Circinus. In this galaxy the point spread function (PSF) was derived 
by comparison to higher resolution VLT/NACO data (see MS06 for details on the fitting 
procedure). A summary of the data, including the specific spatial resolution achieved for 
each galaxy, is given in Table [TJ 



2.3. Extraction of the Emission-Line Gas Properties 

In order to draw general conclusions about the properties of the ionized gas in AGNs, all 
of the galaxies in the sample are evaluated using a consistent method to minimize differences 
due to spatial resolution and the shapes of the spectral profiles. Before discussing our 
findings, we first describe the methods used to fit the spectral line profiles and visualize the 
results. 



2.3.1. Two-dimensional Flux, Velocity, and Dispersion Maps 

We derived emission-line flux distributions, radial velocity and velocity dispersion maps 
of the NLR (narrow Br7 emission) and the CLR using the IDL code LINEFIT. This method, 
described in Davies et al. (2007, hereafter D07), fits the emission lines by convolving a 
Gaussian with a spectrally unresolved template profile (an OH sky emission line) to the 
continuum-subtracted emission line profile at each spatial pixel in the data cube. The con- 
tinuum is estimated using a polynomial function fitted to the line- free regions adjacent to 
the emission line. The continuum and line intervals must be set as input when executing 
the program. To allow for the best possible extraction of the line properties, the width of 
the fitting window is introduced manually after visual inspection of the line profile with 
the visualization tool QFitsVieWLj. The fitting window should cover the complete spectral 
range of the line and enough range of continuum to the left and right of the line profile. 
We estimate that the uncertainty in the continuum subtraction is on the order of 5%. This 
method allowed us to obtain uncertainties for the kinematic maps in the range of (5 — 20) km 
s _1 . The dispersion extracted by this fitting procedure is already corrected for instrumental 
broadening. 

In the central regions of Seyfert 1 galaxies, the Br7 line profile consists of a narrow 



1 see http://www.mpe.mpg.de/~ott/QFitsView/ 
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component (FWHM < 400 km s _1 ) superposed on a broad component (FWHM > 1000 
km s _1 , see e.g. Fig. 1 of D07). In these cases the narrow Br7 properties were extracted 
using the same procedure as described above but utilizing the spectral region of the rather 
flattened top part of the broad component as the interval where the continuum is fitted. 
Since this study is focused on the NLR in Seyfert galaxies, hereafter we refer to the narrow 
component of Br7 just as Br7. 

As it was shown in MS06 for the case of the Circinus galaxy, the [Ca vm] line sits on 
top of the stellar 12 CO(3-l) A2.323/iin absorption bandhead. This occurs in all AGNs with 
[Ca vm] detection, and must be corrected in order to restore the true emission-line profile. 
This correction was performed by convolving stellar templates of red (super)giants with a 
Gaussian broadening function, and varying its parameters to minimize the x 2 ? which was 
measured across the 12 CO(2-0) bandhead (free of emission lines). The 12 CO(3-l) bandhead 
of the template was then convolved with the broadening function and subtracted from the 
galaxy's spectrum (see Fig. 8 of MS06). 

Despite some asymmetries observed in the emission line profiles, we were successful in 
fitting a single Gaussian component to Br7 and the coronal lines (except in NGC 1068, 
see Section 12.3.21) . Typical fits to the [Si vi] emission line profile are shown in Figure [TJ 
The instrumental spectral profile is implicitly taken into account by convolving the assumed 
emission line profile (the Gaussian) with a template line (tracing the effective instrumental 
resolution). As a result, small asymmetries in the emi ssion line profiles ca n be well matched. 



For more details on the key features of LINEFIT see iDavies et al.l (120091 ). This sample rep- 
resents the variety of different profile shapes we find across the IFU fields, and demonstrates 
the high quality of the data and the accuracy of the line-fitting. As can be seen in Figure [TJ 
redshifted or blueshifted tails are common features in the coronal line profiles. An extreme 
special case is observed in the CLR of NGC 4151 (see Figure [2]). While fitting a single 
Gaussian with LINEFIT allows us to reliably measure the [Si vi] properties in the major 
part of the field of view, despite the apparent presence of a broad component in the very 
central pixels which has a considerably lower S/N making its characterization rather uncer- 
tain, the NE and SW regions (regions A and H) exhibit complex profiles in which it is rather 
difficult, if not impossible, to correctly identify the "true" number of components (Crenshaw 
et al. 2000 and Storchi-Bergmann et al. 2010 present studies of the NLR in NGC 4151 
based on two-Gaussian fits to the emission-line profiles). For this reason, and for reasons 
of consistency and simplicity, we decided to limit ourselves to one single main component, 
as in the rest of the galaxies. Furthermore, this approach facilitates the modeling of the 
velocity maps and permits a direct comparison with the fits to the Br7 line profiles, which 
are normally well approximated by a single Gaussian. These fits, in combination with the 
velocity channel maps, which are more sensitive to asymmetries and faint components (see 
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Section r2.3.2p . allowed us to develop a general understanding of both the kinematics and the 
origin of the ionized gas emission. 

In addition, integrated intensity maps of Br7 and coronal lines were constructed by sum- 
ming up all spectral channels containing line emission in each spatial pixel, and subtracting 
from each the continuum level. These measurements result in similar flux distributions as 
those obtained with LINEFIT, confirming that the single Gaussian fits are good approxima- 
tions. 



2.3.2. Velocity channel maps 

The NLR and CLR of NGC 1068 are unique with respect to the other galaxies observed, 
presenting multiple spectral components (see Figure |3]). In consequence, a single- Gaussian 
fit would be insufficient to derive the properties of the ionized gas, and a multi-Gaussian 
approach would be difficult and uncertain due to the large number of degrees of freedom 
to constrain in such fits. In order to take into account the asymmetries (and the multiple 
components as in NGC 1068) in the coronal lines spectral profiles in a consistent manner, a 
velocity tomography of these lines was performed. 

Flux images in different velocity channels were constructed in a similar way to that used 
to create the integrated line flux image (see previous Section). The data cubes were "sliced" 
along the spectral coverage of each emission line into small wavelength ranges to create each 
channel map. The slices were obtained after the subtraction of the continuum determined 
as averages of the fluxes from both sides of the emission line. In most of our data, each 
channel corresponds to ~ 80 km s _1 , approximately the spectral resolution of the data. The 
number of channels depends on the width of the emission line for each galaxy. The zero 
velocity is adopted as the one corresponding to the peak wavelength of H2 1-0S(1) within an 
integrated spectrum centered at the nucleus. In all cases, this corresponds to approximately 
the systemic velocity of the galaxy. 

The velocity channel maps provide different information from that of the 2D velocity 
maps. The latter is dependent on the relative intensities of different kinematic components 
eventually present to form an emission line showing velocities which usually correspond to the 
components with the highest fluxes. In the velocity channel maps, the emission in the wings 
of the profiles is fully mapped, allowing the detection of gas with higher radial velocities, 
but with lower emission intensity. 
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3. General Properties of the Ionized Gas 

This section contains a discussion of the general results of our near-IR observations, 
along with comparisons to relevant measurements at other wavelengths. More detailed anal- 
yses and interpretations can be found in the following Sections. Specific details and analyses 
for individual objects can be found in the Appendix. 



3.1. Spatial Distribution and Kinematics of the Narrow Line Region 

Figures H] to |6] present the Br7 flux, velocity and dispersion maps for the full galaxy 
sample. In all panels, the contours delinate the K— band continuum emission and the position 
of the AGN, defined as the peak of non-stellar continuum emission at 2.2 /j,m, is marked with a 
cross. Regions in white correspond to pixels where the line properties are uncertain and thus 
were masked out. These rejected pixels in the velocity and dispersion maps are those with 
a flux density lower than 5% of the peak of Br7 emission. The Br7 fluxes were measured 
by integrating the flux inside a rectangular region that included all the visible emission, 
F(Br7) int , as well as within a circular aperture equivalent to two times the spatial resolution 
achieved in each galaxy (2xFWHM) and centered at the nucleus, F(Br7) nuc . In NGC 6814 
and NGC 7469 these quantities were measured in both the OSIRIS and SINFONI datasets, 
and the results are consistent within 5%. The measurements are presented in Table |2j 

Despite the diversity of the observed morphologies, most of the Br7 emission consists of 
a marginally resolved core (except in NGC 3783 which is unresolved) and extended emission 
out to several tens of pc. The emission in all cases is diffuse or filamentary, and it is difficult 
to determine whether it further breaks down into compact knots or blobs such as those 
found in Ha or [O ill] images, even though the resolutions are comparable. The only two 
cases where the nuclear emission shows detailed structure (narrow-line clouds) are NGC 1068 
and NGC 4151. In order to quantify the extent of the Br7 emission in all galaxies, we have 
plotted in Figure [7]the azimuthally averaged surface brightness radial profiles of the Br7 flux 
distribution. The azimuthal average of the Br7 emission is consistent with a full-width at 
half-maximum (FWHMb^) less than 60 pc in all of the AGN, with a mean FWHMb I7 for the 
sample of 22 pc (Figure [7] and Table [2]). These measurements reflect the size of the compact 
and bright core in each galaxy. We also measured the extent of the photometric semimajor 
axis of the Br7 emission, Rb^, as well as its PAb^ (Table [2]). The measurements of Rb^ and 
PABr 7 were done directly on the images, using as reference the contours corresponding to 5% 
of the peak of Br7 emission. A comparison of the measured Br7 sizes with measurements 
from [O ill] images is presented in Table |2J The size of Br7 emission in two objects (NGC 3783 
and NGC 6814) is consistent with the sizes of the NLR estimated from [O ill] images. The 



- 12- 



HST [ O ill] image of NGC 7469 s hows a compact bright core (r < 180 pc), also seen in UV 



images (jMunoz-Marin et al.l 120071 ). and extended emission up to r ~ 1300 pc. The OSIRIS 
Br7 flux map of this galaxy shows a similar tendency (see Figure E]). While the compact 
emission is likely associated with the AGN, the extended emission can be attributed to 



circumnuclear starformation (IHeckman et al.lll986l ; iGenzel et al.lll995l ). Therefore, the NLR 



in this galaxy must be very compact and its size is then consistent with the measurements 
of R Br7 (130 < r < 180 pc). In Circinus, NGC 1068, NGC 4151 and NGC 2992, the 
observations are limited by the relatively small FOVs of SINFONI and OSIRIS, but the Br7 
flux maps at these scales show very similar morphologies to those of the NLR as traced by 
[O ill] or Ha images. All these facts confirm the validity of using Br7 as tracer of the NLR. 

As can be seen in the velocity maps of Figures @] to [HI there is a great deal of diversity 
in the kinematics of the NLR at these scales in the sample galaxies. Three cases are clearly 
identified: (1) velocity fields dominated by rotation as in Circinus and NGC 7469, (2) dis- 
turbed rotational patterns (NGC 2992, and NGC 6814), and (3) velocity fields dominated 
by non-circular motions (NGC 3783, NGC 4151 and NGC 1068). As it was discussed in 
Section 12. 3. 2l the Br7 and [Si vi] line profiles of NGC 1068 are composed of multiple spectral 
components. This is clearly affecting the velocity and dispersion maps and therefore they 
cannot be used for analysis. The kinematics of the NLR and CLR in NGC 1068 will be 
analyzed in Section 14.11 by means of the velocity channel maps. 

The velocity dispersion maps, or o maps, show a wide range of values ranging from 30 
km s _1 to ~ 300 km s -1 . Interestingly, in each galaxy, there exist regions where the velocity 
dispersion clearly increases. This will be further discussed in Section 14.11 Table [2] lists the 
mean and the highest a measured for each of the galaxies. 



3.2. Spatial Distribution and Kinematics of the Coronal Line Region 

In general, the S/N of the spectra is high enough to reliably measure the coronal gas flux 
distributions in all galaxies of the sample. As it can be seen in Table [U coronal line emission 
was detected in the seven Seyfert galaxies of the sample. However, the kinematics can only 
be accurately measured where the emission line is strong. Because [Si Vi] is generally the 
most prominent coronal line in the K— band, it was used for the measurements of the size and 
kinematics of the CLR. Figures [8] to [TU] show the [Si vi] flux, velocity and dispersion maps 
for the SINFONI and OSIRIS subsamples. In all panels the contours delinate the K— band 
continuum emission and the position of the AGN is marked with a cross. Once again, in 
the kinematic maps we only display velocities for regions with a flux density higher than 
5% of the peak of [Si vi] emission. The 2D flux, velocity and dispersion maps of the second 
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strongest coronal line in our spectra [Ca vm] look very similar to the equivalent [Si vi] maps, 
so they are not shown here. Measurements of the [Al ix] line also result in similar flux, a 
and velocity distributions. Figure [TT1 demonstrates the similarities between the morphologies 
of [Si vi], [Al ix] and [Ca vm] in the three galaxies where we detected the three coronal 
lines. The quantitative analysis of the line fluxes over the entire CLRs, F([Si Vl])i nt , as 
well as within a circular aperture equivalent to two times the spatial resolution achieved in 
each galaxy (2xFWHM) and centered at the nucleus, F([Si vi]) nuc , is presented in Table |3j 
Once again, in NGC 6814 and NGC 7469 these quantities were measured in the OSIRIS and 
SINFONI datacubes, and the results from both datasets are consistent within 5%. 

As can be seen in Figures [HI to [TUl in all cases the CLR is resolved (except in NGC 3783) 
and the peak of emission is located at the position of the AGN (or very close to it), being 
spatially coincident with the peak of Br7 emission. This supports the hypothesis that most 
of the Br7 in these objects is produced by the AGN, rather than star-formation (probably 
except in Circinus and NGC 7469, see the Appendices IA.1I and IA.7J) . Despite the diversity 



of the observed morphologies, most of the CLRs consist of a bright nucleus and extended 
diffuse emission along a preferred position angle, which usually appears to coincide with that 
of the Br7 gas. In NGC 1068 and NGC 4151 the emission is enhanced at ~ 0.35" north 
and west of the nucleus, respectively. In each of these galaxies, this bright knot is spatially 
coincident with the jet/cloud interaction region (the "northern tongue" of molecular gas in 
NGC 1068, M filler-Sanchez et al. 2009), just as it occurs with the secondary peak of Br7 
(see Section [5.5. II) . 

In order to quantify the extent of the CLR, we have constructed azimuthal averages 
within circular radial annuli of the flux distribution of [Si vi]. The resulting spatial profiles 
are shown in Figure [12j The FWHM of the CLR is less than 60 pc in all of the AGN, with 
a mean FWHMjg; vi ] for the sample of 24 pc (Figure [T21 and Table [3]). These measurements 
reflect the sizes of the compact and bright [Si vi] cores, which are equivalent to the sizes of 
the Br7 cores. We also measured the extent of the photometric semimajor axis of the [Si vi] 
emission (R[si vi])> as we U as hs PArgj vi ], in the same way as it was done for the NLR. The 
results shown in Table [3] indicate that the CLRs and the NLRs are extended along the same 
PA, and that the CLRs are slightly more compact than the NLRs. 

The fact that the flux distributions of [Si vi], [Ca vm] and [Al ix] are resolved (or 
marginally resolved) and extended, support an origin for these coronal lines in the inner 
part of the NLR or in the transition region between the BLR and the NLR, as suggested by 
previous authors (e.g. Prieto et al. 2005). As can be seen in Figure [TT], [Ca vm] is always 
less extended than [Si vi] , and [Al ix] is always the most compact of the three coronal lines. 
Recent studies (Rodriguez- Ardila et al. 2006) support the formation of the observed stratified 
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CLR through photoionization by the central source, with higher ionization potential species 
located closer to the AGN. Our measurements of [Al ix], which has the highest ionization 
potential of the three, are consistent with this hypothesis. However, it is surprising the 
smaller extent of the [Ca Vlll] emission region when compared to that of the [Si vi] , which 
has higher ionization potential. This fact can be attributed to depletion of gas-phase calcium 
as it is deposited in dust grains, which can alter its gas-phase abundance by large factors 
(e.g. Groves et al. 2004). Still, the apparent differences in size might be merely an artifact 
of the lower S/N of the weaker lines. 

As can be seen in the velocity maps of Figures [8] to [TOj, the kinematics of the CLR in the 
majority of the AGNs appears to be dominated by non-circular motions, except for Circinus, 
NGC 2992 and NGC 6814, in which it appears to be dominated by rotation. When comparing 
the velocity fields of [Si vi] and Br7, one notices that the kinematics of the NLR and CLR in 
five of the seven Seyfert galaxies shown in Figures [8] to [10] are very similar, suggesting that 
the mechanism responsible for the gas motions is the same in the two regions. As it will be 
seen in Section [O] this is also true for the NLR and CLR in NGC 1068. The only exception 
is NGC 7469, in which the NLR appears to be rotating and the CLR exhibits a completely 
different velocity field, probably dominated by radial motions. A detailed kinematic analysis 
of the NLR and CLR will be presented in Section HJ 

The velocity dispersion maps shown in Figures [HI to ITUl confirm the fact that the coronal 
lines tend to be broader than the lines in the NLR, but narrower than those in the BLR 
(400 <FWHM< 1000 km s^ 1 ). Surprisingly, the velocity dispersion variations of the gas 
in the NLR are very similar to those observed in the CLR for the majority of the galaxies 
in the sample, providing further evidence for a common mechanism responsible of the gas 
motions in the two regions. Table [3] lists the mean and the highest a measured for each of 
the galaxies. 



4. Kinematic Analysis 

4.1. Evidence for rotation and outflow in the NLR and CLR of AGN 

We have extracted Br7 and [Si Vi] kinematics in the nuclear region of seven Seyfert 
galaxies, allowing for the first time a spatially resolved two-dimensional study of the low- 
and high-ionization gas motions on scales < 300 pc around the AGN. At these scales, the 
Br 7 velocity fields of four galaxies (Circinus, NGC 2992, NGC 6814 and NGC 7469) appear 
rather regular, despite some deviations from pure rotation. A similar situation occurs in 
three velocity maps of [Si vi] emission (Circinus, NGC 2992 and NGC 6814). This fact 
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allows us to use a simple axisymmetric disk model as a first approximation to characterize 
the NLR and CLR kinematics. 

This first chara cterization was done u sing the kinemetry method, developed and de- 



scribed in detail by iKrajnovic et all (120061 ) . Briefly, kinemetry is an extension of surface 
photometry to the higher-order moments of the velocity distribution. Kinemetry was origi- 
nally designed for use with very high S/N (> 100) and relatively extended stellar kinematic 
data. To apply this method to our lower S/N emission line data covering small spatial scales, 
the breadth of the analysis must be somewhat restricted. We therefore employ kinemetry 
in a more limited capacity; rather than using kinemetry to measure and interpret subtle 
kinematic features of a velocity field, we instead use it to determine the strength of devia- 
tions of the observed velocity fields from the ideal rotating disk case. This is identical to 
assuming that any deviations from the ideal case that might occur in a disk (e.g. in- and 
outflows, warps, spiral structure) induce similar power in the higher Fourier coefficients than 
those caused by the noise, and therefore the analysis of these coefficients would be difficult 
and uncertain. We then considered the entire velocity field at once and solved for a global 
position angle (PA) and inclination i. The kinematic center is assumed to be spatially co- 
incident with the AGN defined as the peak in the K— band nonstellar emission. The PA of 
the kinematic major axis is measured counterclockwise from North (0°) to East (90°) and 
have values in the range [—180, 180] degrees. The gradient of velocities, from blueshifts to 
redshifts, determine the sign of the PA. 

The PA and i of the kinematic major axes as listed in Table S] were obtained by min- 
imizing the difference between the SINFONI and OSIRIS velocity fields and the "circular" 
velocity field (the 2D image reconstructed from only the cosine term Bl in the kinemetric 
expansion, Krajnovic et al. 2006). The quality of the fit is represented in each case as the 
average of the absolute values of the velocity residuals map (//), obtained after subtracting 
the rotational component to the Br7 and [Si vi] velocity maps. Residuals greater than the 
estimated kinematic errors (~ 20 km s~ x , see Section f2.3.ip represent either deviations from 
circular motion (bars, inflows/outflows, etc.), or artifacts induced by the noise. For compar- 
ison purposes, we include in Table H] the best fit PA and i of the molecular gas disks derived 
by H09 for the majority of the sample galaxies. These values are in addition consistent with 
those determined from the stellar kinemati c s in t he central < 300 pc and the large scale 



disk, as checked by H09 and iFriedrich et al.l ( 120101 ) for all the galaxies in the sample. Three 



cases are clearly identified in Table H] (note that NGC 1068 is excluded from this analysis). 
Regions with residuals consistent with the estimated kinematic errors (/i < 20 km s _1 ) and 
best-fit PA and % in good agreement with those determined from stellar and/or H 2 kinemat- 
ics (within a typical uncertainty of 12°) are interpreted to be consistent with co-planar disk 
rotation (NLRs in Circinus, NGC 2992 and NGC 7469). The regions with // > 20 km s" 1 
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are likely to be dominated by non-circular motions (NLR and CLR in NGC 3783, NGC 4151 
and NGC 6814; and CLR in NGC 7469). Finally, in two CLRs (Circinus and NGC 2992), 
we note inconsistency between the geometric parameters of the [Si Vl] disk and those of the 
H 2 disk (PA differences greater than the measurement error of ~ 12°, see Table HJ), but the 
residuals indicate a good fit (// < 20 km s _1 ). Since the derived kinematic parameters of 
the CLR in the Circinus galaxy could be affected by the instrumental artifacts present at 
~ 0.25" and ~ —0.25" in the FOV, we consider this result as uncertain. In consequence, a 
detailed modeling of the [Si vi] kinematics in this galaxy was not attempted. The particular 
case of NGC 2992 will be discussed in detail in Section 14.21 

Now we concentrate on the objects with residuals greater than the kinematic errors 
(/i > 20 km s _1 ), indicating the presence of non-circular motions. A second characterization 
was done keeping the best fit parameters of the rotating molecular gas disk constant during 
the kinemetry expansion. However, in this case the measured errors are even larger than in 
the previous characterization indicating that rotation alone cannot explain the kinematics of 
these objects. This can be clearly seen in Figure [TBI were we show the velocity fields of H 2 , 
Br 7 and [Si vi] in NGC 3783, NGC 4151, NGC 6814 and NGC 7469 (the four galaxies with 
// > 20 km s _1 ). The velocity maps of the NLR and CLR in these objects are very similar 
(except in NGC 7469), and present redshifted and blueshifted velocities in the direction 
opposite to the sense of rotation as traced by the H2 velocity maps. In NGC 6814 there is 
a region north of the AGN in which the Br 7 and [Si vi] velocities are redshifted but the H2 
velocities appear approximately at the systemic velocity of the galaxy, or even blueshifted, 
showing a complete different behaviour. As the direction of motion (radial velocity changing 
from blueshift to redshift) occurs along a single PA, there is no evidence of a significant warp 
(e.g. twisting of the kinematic major axis PA and/or inclination angle). Thus, any approach 
to reproduce the observed kinematics by a rotating or warped disk model can be excluded. 
The only exception is NGC 6814 in which there appears to be a twist of the kinematic major 
axis, which could be a signature of either a warped disk or a second kinematical component 
superimposed on co-planar disk rotation. 

In order to better understand such motions, we make use of the velocity channel maps 
described in Section 12.3.21 Figures IT4HT81 show the slices around the [Si Vi] line through the 
data-cubes of NGC 1068, NGC 3783, NGC 4151, NGC 6814, and NGC 7469, respectively. 
For completeness, Figures [19] and [20] show the velocity channel maps of [Si Vi] in NGC 2992 
and Br7 in Circinus, respectively. Since the Circinus datacube is affected by instrumental 
artifacts around the [Si Vi] wavelength, we do not attempt to produce velocity channel 
maps of this line. It should be noticed that in the rest of the galaxies the velocity slices 
of Br7 are very similar to the equivalent [Si vi] maps, so we do not show them here. The 
only exception is NGC 7469, but since the Br 7 velocity map of this galaxy is exceptionally 
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regular (see Figure fl~3]) . the velocity slices do not add significant and useful information to 
this work. 

As can be seen in Figures [T4T - fr8"| the velocity channel maps reveal blueshifted and 
redshifted velocity components on each side of the nucleus, whose speed appears to increase 
with radius to a value that exceeds 200 km s _1 . Note that in NGC 3783 this occurs mostly 
in the north part of the galaxy. These approaching and receding regions are in fact the 
main kinematical components in the radial velocity maps, but the magnitude of the velocity 
vectors is smaller. These radial motions indicate that a considerable amount of ionized gas 
must be moving either towards or outwards the nucleus rather than orbiting it on circular 
paths. The non-circular velocities show values in the range 200-500 km s _1 , except for 
NGC 1068 which exhibits a maximum velocity of ~ 1400 km s^ 1 at ~ 1" from the nucleus. 
The velocities of these bright clouds at ~ lOOpc from the nucleus are too high to be explained 
by any sort of reasonable gravitational potential (including contributions from the SMBH, 
nuclear stellar cluster, and bulge; Das et al. 2007), which rule out infall or orbital motions 
as being dominant. Rotational or gravitational infall models would require unreasonably 
high masses in the central ~ 50 pc (greater than 10 9 M ) to produce the high velocities 
that are seen (> 300 km s _1 ). In the case of NGC 6814, these high velocities suggest that 
a radial component superimposed on disk rotation is responsible for the apparent change 
of the kinematic major axis PA. Furthermore, gas falling towards the nucleus is expected 



to increase its velocity as it approaches the AGN (jMuller-Sanchez et all 120091 ) or present 



defined structure (signatures of bars or inner spirals) in the radial velocity maps or the 
velocity residuals (jSchoenmakers et al.lll997l ; iDavies et al.ll2009[ ). properties that are absent 



in our data. The observed velocity variations as a function of distance provide support for 
a non-gravitational force which is accelerating the gas from the nucleus and therefore can 
only be explained in terms of radial outflow. 

Several other pieces of evidence suggest that the non-circular motions observed in the 
kinematics of the NLR and CLR correspond to outflows (see the Appendix for a discussion 
of the individual galaxies) : 

• For the galaxies with existing radio images (NGC 1068, NGC 4151 and NGC 6814), we 
found that the collimated radio emission is spatially coincident with the NLR/CLR gas 
(or at least the two are aligned, see Figure I2TJ) . The direction of motion seen in [Si vi] 
and Br7 is blueshifted in the jet pointing toward us and redshifted in the counterjet. 
This is a clear indication of a bipolar outflow of material from the nucleus. 

• Interestingly, in most of the galaxies there exist regions where the velocity dispersion 
clearly increases. These increments are observed at the nucleus and/or along the 
directions of the blueshifts and redshifts, providing evidence for radial outflow, as the 
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combination of the line-of-sight (LOS) velocities from clouds located at the front and 
back parts of the outflow increment the width of the spectral profiles. Gas moving on 
circular orbits or falling into the nucleus along eccentric orbits or under the influence of 
a bar or nuclear spiral is expected to show a rather uniform dispersion map (H09), or 
even a low dispersion along it s trajectory towards the nuc leus, since the gas is expected 
to be gravitationally bound ( iMiiller-Sanchez et al.ll2009l ). 



• Some galaxies show evidence for radial acceleration to a projected distance, followed 
by deceleration (NGC 1068, NGC 3783 and NGC 7469). For example, in NGC 1068, 
the maximum blueshift of ~ —1400 km s _1 in Br7 and [Si vi] occurs at ~ 80 pc 
NE from the nucleus, but there are low velocity components at further distances (see 
e.g velocity channel with -1000 km/s), indicating that the NLR/CLR gas accelerates 
outward from the nucleus to a turnover point, and subsequently appears to decelerate. 
Radial acceleration plus deceleration was first noticed in HST /STIS observations of 
NGC 1068 and NGC 4151 (Crenshaw et al. 2000a, 2000b), being interpreted as radial 
outflow. 

The integral-field observations are revealing kinematic signatures of rotation and/or 
outflow in the NLR and CLR of Seyfert galaxies. Therefore, kinematic models incorporating 
these two types of motions were used to reproduce the observed velocity fields. The fitting 
method and the derived model parameters will be discussed in Section 14.21 



4.2. Kinematic modeling 

To quantify the kinematics of the NLR and CLR in the galaxies with significant non- 
circular motions, we have modeled the observed velocities as extra-planar motions with 
biconical geometry, or biconical outflows, superimposed on co-planar disk rotation. The 
biconical geometry was selected because of its adaptability to the present investigation as 
demonstrated by pr evious morphological and kinematical studies of the NLR in S eyfert 
galaxies with HST (ICrenshaw fc Kraemerl |2000| ; ISchmitt et al.ll2003l ; iDas et al.l 120061 ) . Our 
observations also provide support for the use of such geometry. As can be seen in Figure [13] 
and Figures [T41 to [TBI the outflows appear to be bipolar and are following a preferred direction 
with a certain opening angle, characteristics that can be well matched with a bicone. 

The dis k is modeled as a set of rotating coplanar annuli following the cosine law ap- 
proximation (jKrajnovic et al.ll2006l ). The kinematic center is assumed to be coincident with 
the AGN defined as the peak in the nonstellar emission. The disk model has three free 
parameters: the magnitude of the circular velocity v c , the position angle of the major axis 
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of the disk PAdisk, and the inclination of the disk idisk- Fitting the velocity fields with this 
model would be equivalent to use kinemetry to derive the best-fit disk parameters as it was 
done in Section 14.11 

The second component of the mode l is a biconical ou tflow. The modeling of this struc- 



ture is based on the method described in lDas et all (120061 ) . The model consists of two hollow 



cones having interior and exterior walls placed apex to apex at the position of the AGN. 
This point corresponds to the apex of the bicone. The bicone is symmetrical in geometry 
about the apex. Each cone extends from this point to a maximum distance z max measured 
in parsecs along the bicone axis. The radius of the inner and outer walls of the bicone are 
defined by the inner and outer half-opening angles (0; n and 9 out ), respectively. The orienta- 
tion of the bicone in the sky is defined by a pair of angles: PA bicone , and bicone- PA biconc is 
the position angle of the bicone axis in the plane of the sky measured counterclockwise from 
North (0°) to East (90°), and bicone is the inclination of the bicone axis with respect to the 
line-of-sight (LOS), being zero when the bicone axis is perpendicular to the LOS, positive 
when the northern part of the bicone (or eastern part if PAbicone — 90°) is inclined towards 
the observer, and negative otherwise. The velocity at each point is a three-dimensional vec- 
tor moving radially away from the AGN. Since the non-circular (outflow) velocities appear 
to increase with distance r from the nucleus, we tested velocity laws of the form v ex r n , 
and found that a simple proportionality (n=l), or r-law (v a = fcir), provided a good match 
to the data. Some galaxies show evidence for radial acceleration to a projected distance, 
followed by deceleration that approaches the systemic velocity (NGC 1068, NGC 3783 and 
NGC 7469). We therefore incorporated a deceleration term to the velocity law of the form 
Vd = "max — k 2 r, where t> max is the maximum velocity reached by the outflow at a distance 
Tt from the nucleus. The velocities along the bicone walls will be zero at z max . Note that Vd 
does not have any impact on the cases where deceleration is not observed (as in NGC 4151), 
since this term is only used at distances r > r t , and therefore it occurs outside the observed 
emission. In these cases, the exact value of z max is not constrained. 

The bicone model was programmed in IDL as a three-dimensional array filled with ve- 
locities for a given parameter set and velocity law. The velocity value of each element of the 
array corresponds to the component of the velocity vector in the direction of the LOS, as 
would be measured by the observer. An example of such model following an r-law of accel- 
eration and deceleration is shown in Figure [22j The three-dimensional structure and front 
and back projections of two test models with PAbicone = 0° and two inclinations, ibicone = 90° 
and ibicone = 0° are shown in this Figure. For the purposes of illustration, the velocities 
are assigned a color based on the amount of redshift or blueshift. The dark purple color is 
the highest blueshifted (or negative) velocity, and the light red is the highest redshifted (or 
positive) velocity, with intermediate velocities given in gradients. Green represents approxi- 
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mately zero velocity. 

Since the position in the LOS direction of the more luminous clouds in the NLR/CLR 
is not known, for each galaxy we performed three fits: (1) using the LOS velocities in the 
front part of the bicone, (2) using the LOS velocities in the back part of the bicone and 
(3) using the average of the LOS velocities at each pixel. With each of these methods it is 
possible to carry out a two-dimensional kinematic fit of the bicone model plus rotation to the 
velocity maps. In addition, for each galaxy we performed two fits, one using only the bicone 
model, and a second one adjusting simultaneously the disk and bicone parameters. The 
best-fit parameters were found by minimizing the sum of the squared differences between 
the model velocity field and that extracted from the data, taking into account the measured 
errors in each spaxel (x 2 )- The best-fit corresponds to the model presenting the minimum 
X 2 , and the confidence intervals are determined by a range of x 2 values greater than \ 2 



ravni 

defined as Ax 2 = X 2 ~ X 2 m in (ILampton et al.l 119761 1 . Before com puting confidence l evels 



we have renormalized x 2 following a procedure already adopted by iBarth et al.l (120011 ) and 



Marconi et al.l (120061 ). which consists in rescaling the error bars by adding a single constant 
in quadrature to the velocity uncertainties. This constant is chosen for each galaxy such 
that the best fit model has x 2 /dof ~ 1. 

The results of applying this fitting procedure to the four galaxies with significant radial 
motions in their CLR (/i > 20 km s _1 ) are summarized in Table |5j Since the velocity maps 
of the NLR and CLR in these objects are very similar (except in NGC 7469, see Figure IT3~j) . 
the fitting procedure yielded best-fit parameters which are equivalent for the two regions. 
Figures E3H26] show the best-fitting model velocity fields of NGC 3783, NGC 4151, NGC 6814 
and NGC 7469, respectively. In all cases, the best fits were obtained using the average of the 
LOS velocities. In NGC 3783 we adjusted a kinematic model consisting of one single cone 
in the north part of the galaxy plus rotation, as part of the southern cone is hidden by the 
plane of the galaxy (see Section l5T4"j) . For comparison purposes, we present in Tabled the 
minimum x 2 values achieved for three different types of models (without renormalization): 
(1) pure rotation (Xrot); (2) pure biconical outflow (x 2 u t) an d (3) rotation plus biconical 
outflow (Xmt+out)- Although the number of fitted parameters is different in each of these 
three models, the number of degrees of freedom remains practically unchanged, since the 
number of data points in the velocity maps (N) is large. Thus, N was used as the number of 
degrees of freedom in each case (Table |5]). In addition, we show in Figure 1271 the best fit with 
pure rotation for a rotation-dominated case (NGC 6814, x 2 t < Xout)> an d the best fit with 
pure outflow for an ouflow-dominated case (NGC 3783, xiLt < Xrot)- I n the two cases, the 
quality of the fit, represented by the velocity residuals map and the x 2 values (see Tabled]), 
decreases considerably when compared to the best fits incorporating rotation plus outflow. 
As can be seen in Figure |27] specific flow patterns are noticeable in the velocity residuals 
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maps. Particularly interesting are the residuals in an outflow-dominated case, which reveal 
signatures of rotation as traced by the H 2 disk. The NLR/CLR kinematics of NGC 1068 
was modeled using a bicone structure similar to the one described here. However, the fitting 
method was slightly diffe rent, as in this galaxy we aime d to reproduce the velocity channel 



maps and PV diagrams (IStorchi-Bergmann et al.ll2010i ). The modeling of the NGC 1068 



kinematics is described in a separate publication (Miiller- Sanchez et al. in prep.), and here 
we only present the main results (Table [5]), which will be used in the discussion presented in 
Section 

Our kinematic models of biconical outflow plus rotation provide a good match to the 
overall velocity patterns, despite slight local variations in the velocity fields. No obvious 
regular patterns are noticeable in the residuals maps, suggesting that the features observed 
in the residuals are only due to noisy pixels in the data (Figures l2"B"H2T)j) . Furthermore, the 
residuals show mostly values between [—20, 20] km s" 1 consistent with the measured errors 
of the velocity maps. As can be seen in Table [5j in two galaxies the best-fit is obtained 
for models incorporating only biconical outflow (NGC 1068 and NGC 7469) and in other 
two galaxies there are no signatures of deceleration (NGC 4151, NGC 6814). Note that in 
NGC 4151 further observations with larger FOVs are required to conclude on this topic. For 
NGC 1068 the best-fit model parameters are i n agreement with those from previous HST 



spectroscopic studies of the NLR in this galaxy (ICrenshaw fc Kraemerll2000uDas et al.ll2006l ). 
However, in NGC 4151 our estimates of PA bicone , t> max and r t are significantly different to the 
ones derived in previous kinematic studies (Das et al. 2005). We attribute these differences 
to the fact that we are incorporating a rotational component and to the shape of the FOV 
of OSIRIS. Unfortunately, our observations do not fully cover the regions where the outflows 
are observed, so we are probably seeing only the starting and lateral parts of the outflow. 
However, the rest of the model parameters are in good agreement with previous kinematic 
models of the NLR in this galaxy. In all cases where a rotational component is present, the 
PAdisk and idisk values are consistent with those determined from the stellar or H 2 kinematics 



(Table H]) and the large scale galactic disk. This will be further discussed in Section 15.11 

The case of NGC 2992 is now revisited. Since the [Si vi] velocity map of Circinus is 
uncertain (see Section I4TTJ) . this is the only galaxy in which we note inconsistency between 
the geometric parameters of the putative [Si Vi] disk and those of the H 2 disk (PA differences 
greater than the measurement error of ~ 12°, see Table HJ), but \x < 20 km s _1 (Section 14.11 
and Table H]). The orientation of the galactic disk in this galaxy agrees with that of the 
H 2 kinematics (Friedrich et al. 2010), suggesting that although purely co-planar disk ro- 
tation matches reasonably well the [Si vi] velocity field, the CLR kinematics is probably 
more complex. NGC 2992 exhibits blue-wing asymmetric coronal line profiles in the whole 
FOV, reaching values up to ~ —300 km s _1 . This situation is also observed in Circinus and 
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NGC 7469 among the sample galaxies (see Figure [T] and Figure 8 of MS06), and has been 
interpreted as evidence for radial outflow (MS06). Further evidence for outflows comes from 
the increments of velocity dispersion with radius, and from the H 2 kinematics, which are 
explained by incorporating an ouflowing component (Friedrich et al. 2010). Therefore, we 
interpret the variations in the geometrical parameters of the disk as evidence for radial out- 
flow, and decided to use a biconical outflow model superimposed on co-planar disk rotation 
to explain the CLR kinematics in this galaxy. The results of applying this fitting procedure 
to the CLR of NGC 2992 are presented in Table M and Figure [2EJ This model provides also 
a better fit to the Br7 velocity field when compared with the pure rotation model. Further- 
more, the best -fit parameters of the bicone are consistent with previous studies of the NLR 



in this galaxy (jVeilleux et al.ll2001t iGarcia-Lorenzo et al.ll200ll ). 



Table [7] summarizes the main results of the kinematic analysis. In the galaxies where 
the best fit incorporates rotation and radial outflow, the dominant component is considered 
to be the one exhibiting a smaller \ 2 when fitting the velocity field just with that type of 
motion. A visual inspection of the residuals maps was also used to identify the dominant 
component in the NLRs and CLRs. The results of these two approaches are consistent 
with each other. As it was mentioned in Section 13.11 in four galaxies (Circinus, NGC 1068, 
NGC 2992 and NGC 4151) the size of the NLR, as traced by [O ill] emission, is larger than 
our FOVs. However, the model parameters of NGC 1068, NGC 2992 and NGC 4151 are 
consistent with those derived in previous kinematic studies of the NLR in these galaxies 
with larger FOVs. All CLRs are fully covered by the SINFONI and OSIRIS FOVs, except 
in the Circinus galaxy. Since we have not been able to characterize a putative outflowing 
component in Circinus on scales of less than 16 pc from the nucleus, although there exist 
spectroscopic and morphological evidences of its presence (Marconi et al. 1994; Prieto et al. 
2005; MS06), this galaxy is excluded from the discussion in the following sections, except for 
Section 15. 1\ where the presence of rotational components in all galaxies is discussed. 



5. Discussion 

5.1. Inner disks of ionized gas 

In general, the Br7 velocity fields in the sample galaxies exhibit co-planar rotational 
components, being the only exception NGC 1068 (see Table [7]). A similar situation is ob- 
served in the [Si Vl] velocity maps, this time being NGC 1068 and NGC 7469 the only 
exceptions. Furthermore, the best-fit parameters to the NLR and CLR rotational compo- 
nents are consistent with those derived from the stellar/H 2 kinematics. Another important 
fact is that the magnitudes of the rotational velocity vectors in each galaxy are consistent 
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with the values found by H09 for the molecular gas, within ±20 km s _1 . Therefore, the 
dynamical masses in the central r ~ 50 pc estimated from the rotational components of 
ionized gas (considering cr mean of Br7 as the dispersion in the disk, Table 2) are in good 
agreement with those found from H 2 kinematics (between 2 — 20 x 10 7 M ). In addition, 
the consistency of the nuclear kinematics with the large-scale disk orientation reported in 
the literature for the sample galaxies indicates that these inner disks are well aligned with 
the host galaxy. 

Evidence for an underlying rotational component in the NLR is also found in ground- 
based and HST st udies of several oth er objects, such as NGC 1365 (Hjelm & Lindblad 



1996), NGC 2992 flVeilleux et al.l 120011), and NGC 3516 (Arr ibas et al. 1997), and recent 



HST spectroscopic studies of the CLR (JMazzalay et al.ll2010l ) with a general trend for the 



low-ionization gas to be a better tracer of the disk component, while the high- ionization 
gas presents more deviant behavior, usually associated with outflow. Our integral-field data 
confirm this hypothesis and provide accurate measurements of the physical properties of the 
disk and outflow components. 

The presence of inner disks of low- and high-ionization gas provides additional con- 
straints to the classical picture of the NLR in unification models (e.g. Antonucci 1993), 
which postulate that an optically and geometrically thick toroidal structure produces a bi- 
cone of radiation that ionizes the NLR. In this model, the ionized gas is co-spatial with the 
bicone of radiation and is outflowing. Within this context, rotating disks of ionized gas may 
be observed if the ionization bicone intersects the disk plane, causing NLR and CLR emission 
with a rotation s ignature. An exam ple of such interaction has been observed in Mrk 573 



with HST/STIS ( jFischer et al.ll2010l ). This hypothesis can be tested quantitatively with our 
kinematic models. If the angle between the axis of the bicone and that of the host galaxy 
(J3) is smaller than the half openig angle of the bicone ((3 < 6 out ), the ionizing radiation will 
intersect the disk of the host galaxy, and a rotational component will be observed. Because 
the disk has a finite thickness, when /3 > 6 out radial motions are expected to dominate, but a 
rotational component may be present. As can be seen in Table the two galaxies presenting 
only radial outflow in their CLR (NGC 1068 and NGC 7469) have /3 > out , and therefore 
rotational signatures are absent. In NGC 2992, NGC 3783 and NGC 6814 (3 < 9 out , and as a 
result co-planar disk rotation is observed. Finally, in NGC 4151 /3 > 8 out and a rotating disk 
is present, suggesting that either the disk height is high enough to provide a large intersection 
area with the ionization bicone, or there exist hard UV/soft X-ray radiation located outside 
of the ionization bicone in the disk plane. Although the first explanation seems plausible, 
the second should not be ruled out. First, in all four galaxies with rotational and outflow- 
ing components (NGC 2992, NGC 3783, NGC 4151 and NGC 6814; see Table ED the low- 
and high-ionization gas emission extends beyond the borders of the bicone as represented 
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by the models (Figures I23TT261 and Figure [281) . Additionally, in all our sample galaxies, the 
disks of ionized gas are spatially coincident with the disks of molecular gas, suggesting that 
the two types of gas are spatially mixed. The presence of NLR/CLR gas outside the ion- 
ization cones seems to contradict the predictions of the classical torus model of AGN, as 
it would imply the escape of radiation through the walls of the obscuring structure. One 
way to reconcile this discrepancy is to assume that the molecular gas within this nuclear 
region is in a clumpy distribution such that a fraction of ionizing photons escape in different 
random directions outside the bico ne and ionize the gas in the di sk producing the observed 
flux and velocity maps. Recently, iMor. Netzer fc Elitzurl ( 120091 ) analyzed mid-IR spectra 
of AGN using clumpy torus models and found that for a mean number of 5 clouds along 
any equatorial radial line (number of clouds varying between 1 to 10) there is a reasonable 
probability (>~ 1%) to have a LOS that does not encounter any clouds, i.e. UV radiation 
can escape unhindered. The clumpy distribution is also supported by theoretical models 
(jSchartmann et al.ll2005l ; lElitzur fc Shlossmanll2006l ; iNenkova et al.l 120081 ). and observations 
of thermal dust emission (e.g. Jaffe et al. 2004), silicate absorption (e.g. Deo et al. 2007) 
and molecular gas (MS06; H09). 



5.2. Collimation and acceleration of the bipolar outflows 



High degree of collimation ( full opening angle g enerally not more than 15°) is a conspic- 
uous feature of relativistic jets (ILada fc Fichl ll996T ). The ionized gas outflows studied here 
have full opening angles (26 out ) > 50°. In this regime, we consider an outflow collimated 
when its full opening angle is < 70° and less-collimated when it has a value > 90°. 

The outflows observed in the velocity maps of NGC 3783 and NGC 4151 show values 
which increase with distance from the nucleus up to velocities > 200 km s _1 and are colli- 
mated, having half-opening angles (8 out ) in the range 29 — 34° (see Table [5]). The velocity 
channel maps of NGC 1068 (Figure [T41) reveal a similar behaviour, but with a maximum 
velocity which is 4 — 9 times higher than the values estimated in other CLRs and the small- 
est out . These results are consistent with what has been observed previo usly in kinematic 
studies of the NLR of NGC 1068 and NGC 4151 JPas et all I2OO5I . l2006h . In NGC 2992, 
NGC 6814 and NGC 7469 the outflows are less collimated with 9 out between 45 — 57° and 
are not experiencing high accelerations reaching t> max < 200 km s -1 . These results suggest 
the existence of two types of outflows: fast and collimated, and slow and less-collimated. We 
find an anticorrelation between the [Si vi] outflow velocity {v max ) and 9 out of the bicone in 
the sense that more collimated outflows have higher velocities. This relationship is shown in 
Figure [291 Linear fits to the maximum velocities with respect to 6 ont indicate that there is a 
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significant anticorrelation, with Pearson's correlation coefficient of —0.77. The significance 
of this anticorrelation will be investigated further in Section 15.51 

Interestingly, the velocity channel maps of NGC 4151 and NGC 6814 (Figures[T6land [T7|) . 
provide evidence for a combined type of flow: the highest velocities are not only observed 
at the edges of the CLR (or the FOV in the case of NGC 4151), but also at the nucleus. 
These components do not appear in the radial velocity fields because, as it was explained in 
Sect. 12.3.21 the Gaussian fits probe the brightest emission, while the channel maps also probe 
fainter emission (in the wings of the line profiles). This is consistent with the rotational 
component measured in the velocity fields of these two galaxies, as in the vicinity of the 
nucleus the brightest component is the one originating in the galactic disk, while away from 
the AGN the brightest component corresponds to the outflowing one. In the velocity channel 
maps we see high velocity gas at the nucleus, suggesting that the outflow does not leave the 
nucleus at zero velocity. Thus, these observations are probably not revealing acceleration of a 
single wind, but instead the interaction of an already accelera ted wind with the surroun ding 



medium. This scenario i s consi stent with the one proposed by Everett fc Murray! (120071 ) and 



Storchi-Bergmann et al.l (120101 ) who suggested that the interaction of an outflowing nuclear 
wind with the ISM could explain the observed velocities in NGC 4151. A more detailed 
discussion on the origin of the outflows will be presented in Section 15.51 



5.3. Correlation of the outflowing gas and the molecular gas properties 



Unified schemes of AGNs (jAntonuccilll993l ; lUrry fc Padovani 111 9951 ) require an optically 



and geometrically thick structure of gas and dust around the central engine to explain the 
great diversity of AGN classes. The compact sizes (only a few parsecs) determined in recent 
mid-IR interferometric observations (e.g. Jaffe et al. 2004; Beckert et al. 2008; Burtscher 
et al. 2009) require that the obscuring matter be located inside the region where the black 
hole gravity dominates over the galactic bulge. Three fundamental properties characterize 
this structure: (1) anisotropic obscuration so that sources viewed face-on are recognized 
as type 1 objects, those observed edge-on are type 2, (2) dust reemission in the IR of the 
AGN obscured radiation, and (3) collimation power responsible for producing a bicone of 
radiation that ionizes the NLR. In the following we refer to this small structure as torus, 
keeping in mind that its detailed morphology and kinematics are rather uncertain, might not 
resemble closely a rotating smooth torus. As in principle the torus could be the structure 
responsible for the collimation, in this Section we investigate possible correlations between 
the outflowing gas and the molecular gas properties. 

On scales of a few tens of pc from the nucleus, H09 found that the molecular gas in 
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the sample Seyfert galaxies (except for NGC 2992) is in a generally rotating, geometrically 
thick clumpy disk that is associated with the smaller scale obscuring torus. We find an 
anticorrelation between 9 out of the bicone and the molecular gas mass (M gas ) in r < 30 pc in 
the sense that galaxies with greater amount of gas have higher degree of collimation. This 
relationship can be seen in Figure [301 The mass of molecular hydrogen is approximated using 
the typical gas mass fraction of 10% suggested by four independent estimates (H09). The 
anticorrelation coefficient is —0.85. This result confirms that the molecular gas structure seen 
on scales of a few tens of pc has a strong impact on physical processes directly associated 
with the AGN (Davies et al. 2006; Miiller- Sanchez et al. 2009; H09). Furthermore, it 
indicates that the torus is indeed the structure responsible for the collimation of outflows 
on small spatial scales. The accumulation of gas clouds towards the nucleus (higher M gas ) 
increment the collimation power of the torus (stops the outflow expanding sideways). 

A correlation is also found between the [Si vi] outflow velocity f max and the molecular 
gas mass in the sense that galaxies with greater amount of gas have higher outflow velocities. 
This relationship, shown in Figure [3U results from the fact that these two properties are 
anticorrelated with 6 out (Figures [2H] and I3"U|) . The significance of this correlation will be 
investigated further in Section 15.51 



5.4. Orientation effects 

The geometry of the bipolar outflow can be used to test the predictions of the unified 
model of Seyfert galaxies. The basic idea is that the torus outside of the BLR produces a 
bicone of radiation that ionizes the NLR and that Seyfert type depends on viewing angle 
with respect to the bicone axis and half-opening angle of the bicone (ibicone and 9 out ). The 
bipolar outflow in Seyfert 1 galaxies (or intermediate type Seyferts) is expected to satisfy 
the condition | bicone | + #out > 90°, so that a view to the BLR is possible, whereas in Seyfert 
2 galaxies the condition is bicone I + $out < 90°. According to this definition, an AGN with, 
e.g | bicone I = 65°, would be a type 2 nucleus as long as we are not looking into the bicone 
(#out < 25°). Assuming that the plane of the small scale torus is perpendicular to the bicone, 
such an hypothetical AGN would have |ztorus| — 25°, nearly pole-on. We wanted to point out 
this fact to remark the importance of 8 ont in the standard classification of AGN. 

The inclination of the bicone axis in the only Seyfert 2 galaxy studied here, NGC 1068, 
is 9°, nearly perpendicular to the LOS, and 9 ont = 27°. The sum of these two angles is 
< 90° consistent with the prediction of the torus model of AGN for a type 2 object. Since 
0m > Kbiconel, the front part of the bicone is always blueshifted and the back part always 
redshifted, showing blueshifts and redshifts on either side of the AGN. However, blueshifted 
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velocities are expected to slightly dominate the velocity vectors in the NE simply because 
there are more clouds in the LOS moving towards us in this region (the distance across the 
front part of the NE cone is higher than in the back part, due to the inclination). Based 
on the same argument, redshifted clouds should dominate in the SW part of the galaxy. As 
it can be seen in Figure [HJ the bipolar outflow in NGC 1068 is fully consistent with this 
scenario (see e.g velocity channels with —200 and +400 km s _1 ), despite some obscuration 
in the SW produced by the plane of the galaxy (see below). 

Surprisingly, the five intermediate type Seyfert galaxies in our sample (four Sy 1.5 and 
one Sy 1.9, see Table 1) are seen at libiconel+^out < 90° (outside the bicone), which apparently 
contradicts the predictions of the unified model. This puzzling result has also been observed 



previously in kinema tic studies of the NLR of NGC 2992 and NGC 4151 flVeilleux et al. 



200ll ; iDas et al.ll2005l ). Our results, however, are consistent with the unified model within 
the uncertainties. The angle between the outer edge of the bicone and the LOS is defined 
as a = 90° — |ibicone| — ont . Tables [5] and [6] show that in most intermediate type Seyferts 
we have a view to the bicone that is very close to the outer edge, being NGC 3783 the only 
exception. Considering a typical uncertainty of 15° in a (corresponding to the mean of the 
sum of the uncertainties in bicone and 9 out ) a direct view to the BLR is available in most 
cases (except in NGC 3783). It is interesting to point out that other available mesurements 
of a in three intermediate type AGN (NGC 2992: Veilleux et al. 2001, NGC 4051: Fischer 
et al. in prep., and NGC 4151: Das et al. 2005), indicate that the sight line is very close 
to the edge of the cone, rather than "down the middle" (in NGC 4051 the view is actually 
close to pole-on and inside the cone). This is also supported by previous calculations of 
the inclinat ion angles of Seyfert 1 galaxies usin g; results from reverberation mapping and 
X-ray data ( jWu fc Hanll200ll : iBian fc Zhaoll2002l ). Such studies find mean inclination angles 
between 32 — 36°, which translate to bicone = 54 — 58°. Considering the mean half-opening 
angle of the galaxies in our sample 9 out = 40° as a typical value, it is reasonable to conclude 
that most Seyfert Is are actually seen close to the edge of the cone, particularly intermediate 
type Seyfert galaxies (like the ones prese nted in this study, see Table [I]) which are expected 
to have bicone ~ 30 — 45° ( lCurranll2000l ). However, it is intriguing that, according to these 
calculations, many Sy 1.5s would have sight lines that are outside of the bicone, i.e. the BLR 
is not directly seen. Although the uncertainty in our measurements could place most of the 
intermediate type Seyferts in our sample marginally inside the bicone, this may not work in 
the same direction for each one. Further integral-field observations of a larger sample of Sy 
Is and Sy 1.5s are needed to fully test the idea that intermediate type Seyferts are viewed 
at intermediate angles and determine the probability for sight lines in and around the cone. 

Since a is rather small in our galaxies, it seems plausible that, in order to get a direct 
view to the BLR, there exist also small geometric distortions in the torus-bicone system. One 
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possibility would be a small misalignment between the plane of the torus and the bicone axis 
such that they are not completely perpendicular. In this case the condition to observe the 
BLR needs to be modified by adding a correction term: |ibicone| + #out + # > 90°, where 5 
is the angle between ibicone and the true perpendicular to the plane of the torus. Another 
possibility would be that the aperture of the outer part of the torus, where the ionization 
cone escapes, is slightly larger than 6 ont , which in fact is determined by the most inner part 
of the torus. In this case the correction term 5 can be understood as an enlargement to out . 
Even if these two apertures are equal, 9 out can be considered always slightly larger due to 
the existance of a LOS which allows us to see the far side of the BLR. Whatever plausible 
these scenarios may seem, they cannot account for the case of NGC 3783, where our line of 
sight is far outside of the bicone (a = 26°). 

Another way to reconcile the apparent discrepancy between our models and the LOS to 
the BLR (particularly in NGC 3783) is to assume that the collimating structure of molecular 
gas is in a clumpy distribution such that our view to the BLR is almost independent of the 
inclination and opening angles of the bicone. In the clumpy torus scenario, the classification 
of a Seyfert galaxy as type 1 or type 2 depends more on the intrinsi c properties of the torus 



rather than on its inclination. Recently, iRamos- Almeida et al.l (120111 ) fitted models of clumpy 



tori to the spectral energy distributions (SEDs) of seven Seyfert galaxies including NGC 4151, 
NGC 6814 and NGC 7469. Their results do not show a clear trend in the inclination of the 
torus in Sy 1.5 and Sy 2 galaxies (although views of Sy 2s are slightly more inclined than 
those of the Sy 1.5 galaxies), but it is remarkable that for NGC 7469 they found itoms — 83°. 
Their measurements of the inclination angles in the other two galaxies (NGC 4151 and 
NGC 6814) are consistent with ours within the uncertainties. A nalytical models of clumpy 



accret ion disks in combination with SINFONI data were used by IVollmer. Beckert fc Davies 



( 120081 ) to explain the properties and evolution of the torus in six Seyfert galaxies including 
NGC 3783 and NGC 7469. These authors found that while NGC 7469 shows a moderately 
transparent torus, NGC 3783 has the most transparent torus of the galaxies in their sample, 
a result that gives support to our observations. 

The five intermediate type Seyfert galaxies studied here have intermediate inclinations, 
6 U1 <\i\i and | bicone I + #out < 90° (or slightly larger than 90° considering the uncertainties). 
Therefore, they show almost entirely blueshifts on one side and redshifts on the other side 
of the central AGN, simply because there are more clouds with blueshifted (or redshifted) 
velocities at each side of the AGN (one external face and two internal faces of the bicone). The 
intermediate inclination is consistent with unified schemes and may explain the properties 
of intermediate type Seyferts. Recently, Nagao, Taniguchi & Murayama (2000) compared 
several line strengths and line ratios of Sy 1.5s with those of type 1 and 2, and found that 
although the Sy 1.5s have an intermediate value in the line strengths between the Sy Is and 
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Sy 2s, the relative line strengths cover the whole observed ranges of both the Sy Is and the 
Sy 2s. These results suggest that the three types of AGN are basically the same physical 
phenomenon and the differences may be due to inclination effects. In Sy 1.5s a significant 
part of the BLR is obscured by the torus, resulting in composite profiles consisting of both the 
NLR and BLR emission. This is supported by the fact that a spectrum with the appearance 
of an intermediate Seyfert galaxy can be obtained by adding a Seyfert 2 galaxy spectrum to 
a Seyfert 1 galaxy spectrum (Cohen 1983). Alternatively, a clumpy torus could also explain 
the properties of intermediate type Seyfert galaxies. In this case when a cloud is passing 
the line of sight to the central engine, the BLR can be (partially) obscured, resulting in a 
composite profile as well. 

In Figure 1321 we show the geometry of the biconical outflow and the host galaxy disk 
for all galaxies exhibiting radial outflow. For simplicity, the midplane of the galaxy and the 
outer surface of the bicone are portrayed, even though the bicone extends to an inner opening 
angle in the models. The color blue indicates that the majority of the radial velocities in the 
cone are blueshifted, although the back part of the cone in some cases may be redshifted as 
explained above. Accordingly, a red cone indicates that the majority of the radial velocities 
are redshifted. In the case of NGC 1068, the SW side of the galaxy is closer to us, and 
the SW cone therefore experiences more extinction (the red cone). This is consistent with 
the Br7 and [Si vi] images in Figures H] and which shows weaker, and in some places 
absent, emission SW of the location of the AGN. In NGC 2992 the majority of the NW 
cone is blueshifted and the SE cone is redshifted. The SE part of the disk is closer to us and 
therefore the blue cone experiences less extinction, on average, than the red cone. In the case 
of NGC 4151, the SW cone is entirely blueshifted and the NE cone is entirely redshifted. The 
host galaxy disk is close to the plane of the sky, and the SW cone shows less extinction than 
the NE cone shows. This picture is also applicable to the case of NGC 7469 (see Table [5]). 
All the above geometries lead to integrated emission-line profiles that are asymmetric to the 
blue, in agreement with the observed profiles of these galaxies (see e.g. Figure [T]) and the 
majority of previous spectroscopic studies of the CLR based on integrated line profiles. In the 
case of NGC 3783, the PA and inclination of the galactic disk are such that the blueshifted 
cone in the south is in its majority occulted by the disk and the redshifted cone in the north 
experiences less extinction. Thus, this model explains the slightly redshifted tails in the 
emission-line profiles of this galaxy, as well as its radial velocity field. A similar geometry is 
applicable to the case of NGC 6814 (see Table EJ), which also presents asymmetric redshifted 
profiles (see Figures [Tl and ITTj) . These results are consis tent with the an a lysis of the effects 



of extinction by a disk on the NLR profiles presented in lCrenshaw et al.l ( 120101 ) . 
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5.5. Origin of the outflows 

Recently, studies have shown that most of the emission lines seen in the NLR and CLR 
can be accounted for by assuming; a central ionization source illuminating a multicomponen t 
gas (JKraemer et all I2OO0I : broves et aDl2004J JMullaney et al. I bood ; llVfazzalav et al.lboioh . 
Since it is now generally accepted that photoionization is the dominant source of io nization 
of the NLR and C LR gas, with shocks possibly con tributing in localized regions ( jViegas 



20021 ; iPrieto at al.ll2005t iRodrfguez-Ardila et al.ll2006l ). we focus our discussion on dynamical 
properties to gain more insight into the sources of energy that may affect the acceleration of 
the gas. 



5.5.1. Comparison between the gas kinematics and the radio jet 



Figure [21] shows comparisons of the Br7 linemaps to maps at radio wavelengths. The 
outflows are observed in approximate alignment with compact radio emission only in the 
central parsecs of NGC 1068 and NGC 4151, and possibly also in NGC 6814. Note that in 
NGC 6814 the OSIRIS Br7 flux map has three times better spatial resolution than the radio 
image. In the case of NGC 7469, in addition to the diffuse radio emission associated with 
the nuclear starburst r ing, an East- West (E-W) extended nuclear source is also detected 
( lOrienti fc Prietol 120091 ). which, at high spatial resolution, separates into five sources along 
an E-W direction within 0.1" of the nucleus (Lonsdale et al. 2003). Thus, the outflow 
of coronal gas in this galaxy seems to follow the radio distribution at very small spatial 
scales, but a direct comparison between the two data sets could not be made. The 6 cm 
radio emission in NGC 2992 forms a figure-of-eight structure along a PA= 160° (Ulvestad & 
Wilson 1984). In the central 4", the radio continuum is elongated along a PA=38° presenting 
a good correlation with the H 2 and Br7 morphology, but not with the [Si vi] emission, as this 
is very compact (total size of ~ 1.2") and the radio emission present s a cavity in this region . 
In NGC 3783 there exists only an unresolved nuclear radio source (lOrienti fc Prietol 120091 ) 
and in Circinus there is no radio image of comparable resolution available in the literature 
(see Davies et al. 1998 for a comparison of near-IR and radio images in the central 30" x 30" 
of this galaxy) . 

Thus, a perfect spatial correlation between the radio images and the NLR/CLR outflows 
is not seen, only an alignment between them and signatures of interaction between the 
radio jet and the emitting gas in two galaxies: NGC 1068 and NGC 4151. Mundell et 
al. (2003) and Gallimore et al. (2004) suggest that some NLR clouds in NGC 1068 and 
NGC 4151, respectively, are actually responsible for bending the radio jet. Our observations 
are consistent with this interpretation. In these two objects there is a knot of enhanced Br7 



31 



and coronal-line emission at ~ 0.35" north of the nucleus in NGC 1068, and west of the AGN 
in NGC 4151. In each of these galaxies, this bright knot is spatially coincident with one of 
the components of the nuclear radio jet (NGC1068: Component C, Gallimore et al. 2004; 
NGC 4151: Component C, Mundell et al. 2003). The channel maps of Figures [T41 and [TBI 
reveal line flux enhancements at the locations of radio knots with relatively high redshifted 
and blueshifted velocities (speeds in the range 200-600 km s _1 ). These knots can possibly be 
explained by the lateral expansion of hot gas away from the radio jet, as claimed by Axon 
et al. (1998) and Capetti et al. (1997), and indicate that the radio jet has been deflected 
at this location. This flow, however, is restricted to clouds within a few parsecs from the 
nucleus, along the region of the radio jet. The large-scale outflow of the bright clouds cannot 
be accounted for by the radio jet alone and needs a different acceleration mechanism. This 
can be clearly seen in NGC 4151 where the flux distribution at r < 20 pc from the nucleus 
is spatially coincident with the radio jet, but at distances greater than this the radio jet is 
oriented at PA=77° and the bicone at PA=50° (Figure |2T]) . The distinct orientations of the 
two structures suggest that the radio jet is not driving the outflow of ionized gas. A similar 
situation is observed in the rest of the sample galaxies and in NGC 3783 the coronal gas 
clouds seem to accelerate even in the absence of any radio emission. However, at small spatial 
scales, the radio jet may have a contribution to the formation of the outflows, transfering 
momentum to the circumnuclear gas (see Section f5.5.2p . 



5.5.2. An accelerated wind from the AGN 



In the kinematic models proposed for the sample galaxies, the NLR/CLR gas first 
accelerates up to a typical distance of ~ 180 pc (except in NGC 2992, where t> max is reached 
at r = 900 pc) from the nucleus and then decelerates. While the deceleration phase is 
observed in NGC 1068, NGC 3783 and NGC 7469, this stage is not present in the rest of 
our data. 

Dynamical models including radiative forces, gravitational forces and a drag force due 



to th e NLR clouds interacting with a hot ambient medium have been developed (jDas et al. 



20071 ). While these models fail to identify the acceleration mechanism, they do demonstrate 
that drag forces can produce the required deceleration. The possible acceleration mecha- 
nisms, radiative pressure, magnetic fields and thermal winds, fail in explaining the gradual 
acceleration ov er the inner 100 pc, reaching terminal velocities at distanc e s < 10 pc from the 
centr al source (jEverettl 120051 ; IChelouche fc Netzerl l2005t iDas et al.l 120071 : Everett fc Murray 
20071 ). By analogy with the bipolar outflows observed in yo ung stellar objects (Lada & 
Fich 1996; Matzner & McKee 1999), lEverett fc Murray! ((2003) suggested that perhaps the 
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observed war AGN outflows result from the interaction of an already accelerated nuclear 
wind with the ISM, instead of acceleration of a single wind. In protostellar outflows, a 
wide-angled magnetized wind expands like a bubble and sweeps up the ambient medium in 
a momentum conserving manner into a shell at the wind-bubble ambient medium interface. 
The observed velocity la w is a general property of momentu m-conserving; winds blowing 



into 1/r 2 density profiles ( jBally fc Ladalll983l ; IShu et al.lll99ll ). More recently, Matzner & 



McKee (1999) have argued that this velocity relation is recovered for a wide variety of am- 
bient density distributions. Our observations provide evidence for such a two-stage process 
involving the local interaction of high velocity gas (an accelerated wind from the AGN) with 
low velocity gas from the disk (see Section l5\2l) . The AGN winds (e.g. hydromagnetic and 
radiatively driven winds), and probably also the radio jet in some cases, transfer momentum 
to the circumnuclear gas and push the gas from the disk, causing the observed emission-line 
profiles and velocity structure along the bicone. 

In this scenario, the collimation power of the torus and the relative angle between the 
outflow axis and the disk of the galaxy (/3), are crucial to explain the diverse NLR/CLR 
kinematics observed. The correlation between molecular gas mass and v max (Figure I3TJ) indi- 
cates that the collimation power of the torus influences considerably the physical properties 
of the outflows. When the collimation power of the torus is high (high M gas ) and (3 > 6> out , we 
observe fast outflows with no rotational component as in NGC 1068. When the collimation 
power is high and (3 < 8 out we see relatively fast outflows and rotation as in NGC 3783. In 
the cases where the collimation power is low and (3 < 9 out , rotation and outflows with low 
accelerations are present as in NGC 2992 and NGC 6814, and when the collimation power is 
low and (3 > 6 out , we only see a slow outflow with no rotational component as in NGC 7469. 
In addition, as it was discussed in Section [57T| a small fraction of the NLR/CLR gas rotating 
in the disks may be attributed to ionizing photons which have travelled through a clumpy 
distribution of molecular gas. This is probably the case of NGC 4151, in which fast outflows 
(high collimation power) and rotation are observed but (3 > 9 out . 

D07 showed that there is considerable emission from young stars (with ages in the range 
10-300 Myr) within the central few parsecs of an AGN-hosting galaxy. It is thus reasonable 
to consider that the kinematics of the NLR/CLR can be affected by these (post-)starbursts 
at radii approaching a few tens of parsecs from the central engine. Models of starburst- 
driven galactic outflows whose dynamics depends on both radiation and thermal pressure 
from supernova (SN) explosions have been succesf ul in explaining la rge-scale outflows of 



ionized gas in starburst and ultraluminous galaxies (INath fc Silkil2009l ). Using models that 
synthesize the evolution of populations of massive stars we can estimate the mass and energy 
returned from starbursts (e.g. Leitherer et al. 1999, Starburst99). In the solar-metallicity 
case where the mass-loss rate and mechanical luminosity are constant beyond ~ 40 Myr 
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and scale with SFR, E = 7.0 x 10 41 (SFR/M Q yr^ 1 ) erg s" 1 JVeilleux et al.ll2005h . We note 
that for instantaneous star formation the mechanical luminosity is approximately constant 
with a value ~ 10 40 erg s _1 . For a typical SFR of 1 M yr _1 in the central 200 pc of these 
galaxies (D07), the estimated injection of energy from the post-starbursts is in the range 
10 40 — 7 x 10 41 erg s -1 , which is on the order of the estimated kinetic power of the observed 
outflows in four galaxies (see Table [8]), indicating that they are in the limit of kinetic power 
to be able to drive the outflows. Note that in actively starforming galaxies the measured 
SFRs are > 5 Mgvr- 1 . 

Two considerations argue against starburst- driven outflows: timescale and morphology. 
The dynamical timescales of the outflows (tdyn = r t /v max ) range from 0.1 to 4 Myr (see 
Tables El and [H]) • Even considering a longer active phase of the outflow, tipically 10 Myr, 
these timescales are inconsistent with the derived starburst ages but in agreement with the 
typical timescale of the active phase of an AGN. Regarding their morphology, starburst- 
driven outflows observed in nearby starburst galaxies (M82: Westmoquette et al. 2009, 
NGC 253: Muller-Sanchez et al. in prep.) have a different geometry to the one observed 
in Seyfert galaxies, which consists basically of truncated hollow cones. This reflects the fact 
that they originate in the central hundred pc of the starburst galaxy. As the geometry of the 
outflows in our sample Seyfert galaxies is a bicone having a sharp apex, the collimation is 
occuring on very small spatial scales, supporting an AGN origin. We therefore conclude that 
if the central post-starbursts have a contribution on the acceleration of ionized gas in Seyfert 
galaxies, these contributions are modest, being the AGN the major driving mechanism. 



5.6. Mass outflow rates and energetics 

Recent simulations of mergers and quasar evolution have invoked AGN induced outflows 
to explain the correlations between the masses of black holes and the galaxy bulge properties 
(e.g. Silk & Rees 1998; di Matteo et al. 2005). The majority of these models require that a 
significant fraction of the radiated luminosity L^ \ of the AGN couples ther modynamically to 



the su rrounding gas (~ 5%, di Matteo et al. 2005), although recent work by lHopkins fc Elvis 



(120101 ) demonstrate that a two-stage feedback model, similar to the scenario proposed in 
Section I5.5[ requires a much lower energy injection (~ 0.5%Lt, o i) All current simulations 
have been developed for galaxies in which the black hole (BH) is growing rapidly via merger- 
induced accretion, i.e. the feedback is driven by quasars or radio galaxies. Seyfert galaxies 
are generally less luminous tha n quasars and they are likely to be fuelled via secular processes 



associated with disk evolution fjDavies et al.ll2010l ). Although the mechanism responsible for 



triggering accretion on to massive BHs in Seyfert galaxies and quasars may be different, the 
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two lead to the formation of an AGN. The subsequent physical processes associated to the 
AGN phenomenon are expected to be similar in local (Seyfert) and more distant (quasars) 
galaxies. Therefore, it is important to investigate whether the properties of the ionized gas 
outflows in these local AGNs are consistent with the assumptions of the quasar feedback 
models. 

The feedback from the AGN - in the form of the mass outflow rate (M out ) and its kinetic 
power (E out ) - can be estimated using the ionized gas properties and geometry of the outflow. 
We concentrate on CLR gas, in contrast to previous studies in which these quantities were 
estimated using NLR lines (e.g. Holt et al. 2006; Storchi-Bergmann et al. 2010), because the 
outflow properties can be better constrained in the kinematics of the high-ionization species 
(see Figure d3] and Section H~2"j) . but remark that the analysis also applies for the NLR gas 
as the geometry and kinematics of the outflow are very similar in both regions, except for 
NGC 7469 in which only the CLR is in an outflow (see Table [7]). 



For a steady-state outflow, the mass outflow rate is given by 

M out = 2m p N e w max A f 



(1) 



where m p is the proton mass, N e is the ionized gas density, A is the lateral surface area 
of the outflowing region (one cone) and / is the filling factor. The factor 2 is introduced 
to take into account the bicone geometry. We assumed typical values of N e = 5000 cm -3 
and / = 0.001, based on previous studies of the CLR which indicate 10 3 < N P , < 10 4 cm" 3 
and f < 0.01 JMoorwood et all Il99fit loiival Il997l : iRodriguez-Ardila et all boOfit [Schneider 
20061 ). It is important to point out that since the critical density of the CLR is 2 — 3 orders 
of magnitude higher than that of the typical NLR gas (Fergusson et al. 1997; Rodr iguez- 
Ardilla et al. 2011), N e could be up to 200 times higher. However, since N e oc f~ 1 ^ 2 ( Oliva 
19971 ). the filling factor would be very much reduced in those cases. If this calculation was 
performed for t he NLR, N K would be smaller 10 2 < N K < 10 3 cm" 3 and f would be la rger 



0.01 < / < 0.1 flTavlor et al. 112003 : Bennert et al. 1120061 : IStorchi-Bergmann et all 120101 ). A 
and f max are obtained directly from our kinematic models, both measured at the position 
of r t . Therefore, we are confident that the uncertainties in these two last parameters (the 
measurements errors) are negligible compared to those in the assumed values of gas density 
and filling factor. 

The results for the sample galaxies are presented in Table [HJ Our calculations yield mass 
outflow rates ranging from 1 to 10 M yr _1 , except for NGC 2992 which has M out rj 120 M Q 
yr^ 1 . Due to uncertainties in the assumed parameters (N e and /) the derived mass outflow 
rates are probably accurate to no better than about ±0.5 dex. The high M out in NGC 2992 
reflects the fact that this galaxy has the largest lateral surface area of the outflow. However, it 
is unlikely that the radius of the outflow region has been overestimated (r t = 900 pc) since it is 
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consi stent with the one of the figure-of-8 observed in radio observations (jUlvestad &: Wilson 
19841 ) and shorter than the one measured in morphological and spectr oscopic studies of 
Ha and [O ill] emission (jVeilleux et al.ll2001t iGarcfa-Lorenzo et al.ll200l[ ). If anything, the 
outflowing region may be slightly larger. Therefore, the outflow in this galaxy presents 
different characteristics to the ones observed in the rest of the sample galaxies, probably 
having an hybrid origin (stellar and the AGN, see Appendix IA.3J) . In fact, the estimated 
M out in this galaxy is perfectly consistent with those detected in powerful starburst sources 
such as ULIRGs (10 to 1000 M yr" 1 , Veilleux et al. 2005). 



These values can be compared with those reported in the literature. I Crenshaw fc Kraemer 

(120071 ) obtain M out w 0.16 M Q yr -1 fo r NGC 4151 using blueshifted absorption lines in the 
UV spectrum of this galaxy. However, IStorchi-Bergmann et al.l (120 lOl ) report M out w 2.4 M 



yj.-i £ or nqq 4151 based on AO-assisted near-IR integ ral-field observations obtained with 
the instrument NIFS at the Gemini North Telescope. IVeilleux et al.l (120051 ) report values 
of M out calcula ted from warm ionize d gas masses in the NLR in the range 0.1 to 10 M 



yr l . However, iBarbosa et al.l (120091 ) obtained mass outflow rates 10-100 times smaller for 



a sample of Seyfert galaxies observed with seeing-limited integral field spectroscopy. 

We can also compare the measured outflow rates with the mass accretion rates M acc 
necessary to feed the central massive black holes in the galaxies of the sample. The mass 
accretion rate at scales down to one Schwarzschild radius Rs can be estimated from the mass- 
to-luminosity conversion efficiency of a black hole Lboi = f]c 2 M acc , where r\ is the accretion 
efficiency and L^ Q \ is the bolometric luminosity of the AGN. A typical accretion efficienc y 
value for a standard geometrically thin accretion disk is r) = 0.1 (jShakura fc Sunyaevlll973l ). 
The values of Lboi and the resulting mass accretion rates are listed in Table We can 
conclude that, on average, there is ~ 10 2 — 10 3 times more mass outflowing from the nuclear 
region than being accreted into the black hole. A larger outflow rate can be understood if the 
outflowing gas does not originate in the nucleus, but is actually gas from the circumnuclear 
ISM being pu s hed a way by a nuclear wind, consi s tent w ith the two-stage feedback model. 
Veilleux et al.l ( 120051 ) and IStorchi-Bergmann et al.l ( 120101 ) report similar ratios between mass 
outflow rates and nuclear accretion rates in AGN. 

In order to estimate the impact of the ionized gas outflow on the circumnuclear gas in 
the bulge of the host galaxy, it is also important to estimate its kinetic power and compare 
it with the accretion power (Lboi)- Assuming that the relatively large linewdith of the 
outflowing gas reflects a turbulent motion that is present at all locations in the outflow 
region, the kinetic power is E out = l/2M out (w r 2 nax + a 2 ). The resulting values are also 
listed in Table EJ The ratios between the estimated kinetic power and Lboi are in the range 
10~ 4 < Lo U t/Lboi < 5 x 10~ 2 . We note that E ont in half of the sample galaxies exhibiting 
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radial outflow (NGC 1068, NGC 2992 and NGC 4151) fulfill the requirements of the two- 
stage feedback models to be thermally coupled to the circumnuclear gas (E ont ~ 0.005L bol ). 
Furthermore, the fact that none of the measured parameters appear to systematically depend 
on the AGN luminosity (and therefore on the accretion rate), supports the hypothesis that 
the outflow is not a single wind originating in the AGN, but is actually a product of the 
interaction of an already accelerated wind with its ambient medium. It is also likely that the 
kinetic energy of the outflow in these three objects blows away the molecular gas reservoir in 
the central 500 pc necessary to trigger AGN and starformation activity. Assuming that the 
typical lifetime of the outflow is 10 7 yr, the total kinetic energy injected by the outflows to 
the ISM ranges between 3 x 10 56 to 10 57 erg. Assuming that these galaxies have a total mass 
Mtotai ~ 10 10 M Q and a typical gas mass M gas ~5x 10 8 contained within a radius of 500 pc 
(Schinnerer et al. 2000; Garcia-Burillo et al. 2005; Storchi-Bergmann et al. 2010; Friedrich 
et al. 2010), the gravitational binding energy of the gas is E hind ~ CM tota iM gas /r k5x 10 56 
erg. Thus, the outflows may remove a considerable amount of warm/cool gas from the central 
regions of the host galaxy. 

The estimated kinetic power in the three remaining galaxies exhibiting outflows of coro- 
nal gas (NGC 3783, NGC 6814 and NGC 7469) is at least one order of magnitude less than 
might be expected on the basis of the quasar feedback models, even considering a two-stage 
feedback model. Furthermore, these outflows are also incapable of removing a significant 
fraction of the warm/cool gas from the central region since their kinetic energy is less than 
3 x 10 55 erg. As can be seen in Table El and following the discussion presented in Sec- 
tion 15.21 our results suggest the existence of two types of outflows: fast and collimated with 
a relatively high E out , and slow and less collimated with a much smaller E ont - Note that the 
outflow in NGC 2992 does not fit in these scenarios, since it presents characteristics similar to 
those observed in ULIRGs as discussed above. The two galaxies with the highest E out /L ho i 
(NGC 1068 and NGC 4151) are also the only objects in which a clear interaction of the 
surrounding medium with the radio jet is observed. This may suggest that Seyfert galaxies 
with strong and collimated radio emission are also hosts of powerful outflows of ionized gas. 
Further radio and near-IR integral-field observations with high spatial resoultion would be 
required to confirm this hypothesis and to gauge the real impact of ionized gas outflows on 
the evolution of Seyfert galaxies. 



6. Summary and Conclusions 

SINFONI and OSIRIS near-IR AO-assisted integral field spectroscopy has been used in 
this work to probe the kinematics of the NLR and CLR in nearby AGNs. Thanks to our high 
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S/N integral field data we were able to extract the 2D distribution and kinematics of Br7 
and the high-ionization line [Si Vl] in the central < 300 pc with spatial resolutions ranging 
from 4 to 36 pc. Our main results and conclusions can be summarized as follows: 

• In all cases the NLR and CLR are resolved (except in NGC 3783), and their mor- 
phologies are similar, indicating that they are being produced by the same ionizing 
source. The two regions present a bright compact core and extended emission. While 
the size of the NLR seems to be larger than our FOV in some objects, the CLR is more 
compact and could be well measured in the sample galaxies (except in Circinus) with 
sizes ranging from r = 80 pc in NGC 6814 to r = 150 pc in NGC 1068. 

• The kinematics of the NLR at these scales can be put into three groups: (a) velocity 
fields dominated by rotation, (b) disturbed rotational patterns, and (c) velocity fields 
dominated by non-circular motions. For the galaxies with velocity fields dominated by 
rotation, the NLR and stellar velocity fields are similar and indicate that the majority 
of the NLR gas is distributed in a disk. 

• The velocity and dispersion maps of the NLR and CLR in the sample galaxies are 
similar, with a trend for the low-ionization gas to be a better tracer of the disk compo- 
nent, while the high-ionization gas presents more deviant behavior with velocity fields 
dominated by non-circular motions. 

• Several pieces of evidence suggest that the non-circular motions observed in the velocity 
fields of the NLR and CLR correspond to outflows. The main argument favouring this 
interpretation is the evidence for radial acceleration to a projected distance of ~ 100 
pc (followed by deceleration in some cases) observed in the velocity channel maps and 
the radial velocity maps. Other evidences are: (1) For the galaxies with existing radio 
images, the velocity channel maps reveal blueshifts in the direction of the radio jet 
pointing towards us and redshifts in the direction of the jet pointing away from us, 
(2) the velocity dispersion clearly increases at the nucleus and/or along the directions 
of the blueshifts and redshifts and (3) the non-circular velocities are too high to be 
explained by any sort of reasonable gravitational potential. 

• Kinematic models consisting of rotation and/or biconical outflow have been successfully 
applied to the data. The models provide a good match to the overall velocity patterns, 
despite slight local variations in the velocity fields. 

• Regarding the rotational component, the data clearly shows that the disks of low- and 
high-ionized gas are spatially coincident with the inner disks of molecular gas, suggest- 
ing that the three types of gas are spatially mixed. Since high excitation emission is 
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observed outside the borders of the bicone, the torus must be a clumpy medium such 
that a fraction of ionizing photons escape in different directions and ionize the gas in 
the disk. 

The maximum velocity of the outflowing component ranges between 120 km s _1 in 
NGC 6814 to 1900 km s" 1 in NGC 1068, and it is reached at a typical distance of 
~ 180 pc from the AGN. We observe two types of outflows based on their opening 
angle: collimated with full opening angles < 70° and less collimated with full opening 
angles > 90°. 

We find an anticorrelation between the outflow velocity (i> max ) and the half-opening 
angle of the bicone (6 on t), and another anticorrelation between out and the molecular 
gas mass (M gas ) in r < 30 pc. These anticorrelations suggest that galaxies with greater 
amount of gas in their centers have higher degree of collimation (small 8 ont ) and higher 
outflow velocities. This indicates that the torus invoked by unified models is indeed 
the structure responsible for the collimation of the outflows on small spatial scales 
and that the accumulation of gas clouds towards the nucleus (higher M gas ) increments 
its collimation power and produces outflows with higher velocities (stops the outflow 
expanding sideways). 

The observed bipolar outflows are consistent with the torus model of AGN (except in 
NGC 3783). The five intermediate type Seyfert galaxies studied here (four Sy 1.5s and 
one Sy 1.9) are dominated by blueshifts on one side and redshifts on the other side of 
the AGN, consistent with an intermediate viewing angle, whereas the Seyfert 2 galaxy 
NGC 1068 shows blueshifts and redshifts on either side of the AGN, consistent with 
the unified model prediction of a bicone axis near the plane of the sky. Considering the 
uncertainties in the kinematic modeling, and possible misalignments in the torus-bicone 
system, it is likely that we have a direct view to the BLR in most of the intermediate 
type Seyferts (except in NGC 3783). We conclude that the majority of these objects 
(and possibly most Seyfert 1 galaxies) are actually seen close to the edge of the cone, 
rather than down the middle. Alternatively, a clumpy torus may explain the presence 
of broad lines even if the observer is not looking into the bicone, which is likely the 
case in NGC 3783. 

We have estimated mass outflow rates in the range 1 to 10 M yr _1 , which are ~ 10 2 — 
10 3 times the accretion rate necessary to feed the AGN. The estimated kinetic powers 
of the outflows are 10 2 — 10 4 times smaller than the AGN bolometric luminosities. 
We note that E out in half of the galaxies presenting outflows (NGC 1068, NGC 2992 
and NGC 4151) fulfill the requirements of two-stage feedback models to be thermally 
coupled to the circumnuclear gas (E ont ~ 0.005L bol ), and match the M BH — a* relation. 
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The two galaxies with the highest E out /L^, \ (NGC 1068 and NGC 4151) are also the 
only objects in which a clear interaction of the surrounding medium with the radio jet 
is observed. This may suggest that Seyfert galaxies with strong and collimated radio 
emission are also hosts of powerful outflows of ionized gas. 

The high velocity components detected at the nucleus of some galaxies and the high 
M out /M acc ratios indicate that the outflowing gas is not likely a single wind from 
the AGN, but the result of a two-stage process in which the circumnuclear ISM has 
been pushed away by an already accelerated AGN wind. The radio jet may also have 
a contribution to the formation of the outflows at small spatial scales, transfering 
momentum to the circumnuclear gas. 
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A. Description of NLR and CLR gas in individual objects 

In this Appendix, we discuss the emission-line properties of each object, and give a 
brief review of some data in the literature. For each galaxy, we discuss the structure and 
kinematics of the NLR and CLR, and point out relations of these structures with the radio 
jet and [O ill] emission. 



A.l. Circinus 

Circinus displays several characteristics of a typical Seyfert 2 nucleus: the observed line 
ratios (Oliva et al. 1994), an ionization cone observed in [O ill] (Marconi et al. 1994) with 
the corresponding countercone appearing in the [Si Vll] line (Prieto et al. 2005), narrow 
prominent coronal lines in its optical/near-IR spectrum (Oliva et al. 1994; Maiolino et al. 
2000), and broad (FWHM> 3300 km s _1 ) Ha emission detected in polarized light (Oliva et 
al. 1998). 

Br7 and the coronal lines present a relatively symmetric flux distribution extending out 
to a radius of 8 pc along a PA of 0°. Note that the complete extension of the NLR and CLR 
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is not accesible in our observations due to the limited FOV of SINFONI. At larger scales, 
Circinus shows low- and high-ionized gas extending along the minor axis of the galaxy out 
to r ~ 30 pc, with a morphology that is reminiscent of the ionization cones seen in other 
Seyfert galaxies (Marconi et al. 1994, Prieto et al. 2005). 

The Br7 kinematics in the central 16 pc are consistent with ordered rotation in a disk 
with a PA of 18° and an inclination of 65°. Within the errors, these parameters agree with 
the large-scale galactic disk (Freeman et al. 1977). As already pointed out by MS06, this 
suggests that a considerable amount of Br7 in this galaxy must have a stellar origin. This 
is supported by the fact that there exist offsets between the photometric centers at different 
wavelengths and between the peaks of several emission lines (Marconi et al. 1994; MS06), 
indicative of the presence of central super star clusters (MS06; Miiller-Sanchez et al. 2010). 
The velocity field of [Si vi] is also relatively regular, with a best-fit PA of 0° and inclination 
of 57°. However, the best-fit PA to the [Si vi] velocity map is not consistent with that of H2 
or Br7. Furthermore, Circinus exhibits asymmetric and broadened spectral profiles, which 
can be well fitted by a strong narrow component and a weaker blueshifted broad component. 
These two facts are indicative of outflows. As the velocity field is affected by artifacts at 
~ 0.25 and —0.25" caused probably by illumination effects on the SINFONI detector, we 
decided not to perform a kinematic modeling of the CLR in this galaxy. A study of the CLR 
kinematics in Circinus, based on integrated spectra, is presented in MS06. 



A.2. NGC 1068 

NGC 1068 has the most extensively studied NLR in the literature. It was the first one 
to be shown to have a conically shaped NLR (Pogge 1988a). The HST [O ill] images of this 
galaxy were discussed in detail by Evans et al. (1991), Capetti et al. (1997) and Bruhweiler 
et al. (2001). These images show that most of the emission originates in a V-shaped region 
with opening angle of 50°, extending for approximately 750 pc along PA 35°, and 430 pc in 
the direction perpendicular to the bicone axis. The direction of this emission is coincident 
with the radio jet observed by Wilson & Ulvestad (1982). The radio emission is extended 
along PA ~ 33° on the larger scales (Wilson & Ulvestad 1982) but is aligned along the N-S 
direction in the nuclear region (Gallimore et al. 1996). The Br7 flux map shown in Figure H] 
has the same morphology of the [O ill] emission in the central 300 pc. 

The coronal lines in this galaxy have also been extensively studied. Ground-based long- 
slit spectroscopy in the optical and near-IR revealed complex spectral shapes, dominated 
by extreme broadening in all high-ionization lines (Rodriguez- Ardila et al. 2006). The 
HST/STIS spectra show several coronal lines in the optical, from [Ne v] and [Fe Vlll] lines, 
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whose ionisation potential are ~ 100 eV, to the extr eme ionisation [S xn] line, with an 



ionisation potential of 504.7 eV (JMazzalay et al.ll2010l ). This very high-ionization line has 



been reported only in NGC 1068 (Kraemer & Crenshaw 2000). 

The strongest coronal lines in our sample are observed in NGC 1068. The three coronal 
lines contained in the spectral range of our observations, [Si vi], [Al ix] and [Ca vm], 
dominate the spectra in the central region. In our observations, the CLR tends to follow the 
same flux distribution as Br7 and [O ill], but is slightly more compact. Both, the NLR and 
CLR emission are stronger towards the NE suggesting strong extinction towards the SW. 
The kinematics of the two regions are dominated by non-circular motions that can only be 
explained in terms of radial outflow. The evidences are the following: (1) the velocity channel 
maps reveal blueshifts in the direction of the radio jet pointing towards us and redshifts in 
the direction of the jet pointing away from us, (2) there are clear signatures of acceleration 
along a PA of ~ 30°, and (3) the velocities are too high to be explained just by gravitational 
forces. 

The maximum velocity of the outflow t> max is reached at a projected distance of ~ 80 
pc. After this point, the gas appears to decelerate. Due to the limited FOV of SINFONI, we 
do not see the complete deceleration phase, but it has been observed that it reaches systemic 
velocity at z max = 400 pc (Das et al. 2006). The method used to model the kinematics 
is presented in Miiller- Sanchez et al. (in prep.). The kinematics can be well modeled with 
radial outflow along a biconical geometry. The parameters of the bicone are listed in Table [5J 
Furtermore, the kinematics are consistent with the predictions of the Unified Model for a 
type 2 object. The inclination of the bicone axis is 9°, nearly perpendicular to the LOS. 
Based on this value and the external and internal half-opening angles of the cone (0 O ut — 27°, 
9 in = 14°), blueshifted velocities are expected to dominate the velocity vectors in the NE 
simply because there are more clouds moving towards us in this region. Based on the same 
argument, redshifted clouds should dominate in the SW part of the galaxy. In addition, 
redshifted velocities should also be present in the NE corresponding to the back part of 
the outflow, and blueshifted velocities from the front part of the cone should appear in the 
SW. These properties are evident in the velocity tomography of NGC 1068 (Figure [14"]) . 
This galaxy exhibits the fastest and most collimated outflow of the galaxies in our sample 
(E out /L ho \ = 0.05) indicating that the energy released by the outflows can be thermally 
coupled with the ISM. 
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A.3. NGC 2992 



The host galaxy disk is highly inclined at i ~ 65°, which makes it suitable for extra- 
planar emission studies. The galaxy has a companion NGC 2993, which is located 3' SE, 
connected with a tidal bridge, and whose tidal forces might have an important effect on 
the kinematics. The 6 cm radio emission is concentrated in the central parts, forming a 
figure-of-eight structure at PA= 160°. (Ulvestad & Wilson 1984). The [O ill] emission is 
aligned with the galactic minor axis (PA= 120°) and approximately with the radio struc- 
ture. Multi-component optical emission lines were detected, and the velocity field of Ha was 
suggested to be a superposition of rotation and an outflow (Veilleux et al. 2001). 

This galaxy has a compact CLR (r = 100 pc) but an extended NLR (r > 280) pc. The 
two are elongated along a PA~ 30°. It presents several similarities with the Circinus galaxy 
which include: the inclination of the galactic disk, a hard nuclear continuum which produces 
[Si Vl], [Al ix] and [Ca Vlll] emission, the shape of the coronal lines profiles (asymmetric 
with blue wings, see Figure [TJ, extended [O ill] and Ha emission along a PA~ 120° with a 
morphology that is reminiscent of ionization cones, inconsistency between the best-fit PA of 
the [Si vi] velocity field (when fitted with kinemetry) and the PA of the H2 inner disk, and 
the presence of an outflow of molecular gas (MS06, Friedrich et al. 2010). All these facts sug- 
gest that although the [Si vi] velocity field appears rather regular, an outflowing component 
must be part of the kinematics of the NLR and CLR. The Br7 and [Si vi] velocity fields are 
best fitted with a model incorporating radial outflow plus rotation. The resulting parameters 
of the model are presented in Table |6j Interestingly, in this galaxy the outflow reaches its 
maximum velocity of 200 km s _1 at r = 900 pc, the largest distance among the galaxies 
in our sample. Furthermore, it also has the largest opening angle (0 out = 57°). The radius 
of the accelerating ga s in this galaxy is consist ent with the one of the figure-of-8 observed 



in radio observations (jUlvestad fc Wilson! 1 19841 ) and slightly shorter th an the one measured 



in morphological and spect roscopic studies of Ha and [O ill] emission ( 1 Veilleux et al .112001 



Garcfa-Lorenzo et al.l 120011 ) . Therefore, the outflow in this galaxy presents different charac- 
teristics to the ones observed in the rest of the sample galaxies. We interpret these properties 
as evidence for an hybrid origin of the outflow (stellar and the AGN). This interpretation 
is supported by the fact that NGC 2992 is a combination between an intermediate Seyfert 
and a starburst galaxy (Glass 1997). Furthermore, the outtflow in NGC 2992 is similar to 
those observed in nearby starburst galaxies, which usually extend to radial distances of kpc 
(M82: Westmoquette et al. 2009, NGC 253: Miiller- Sanchez et al. in prep.). Aditionally, 
the large 8 out may be an indication of a truncated cone instead of a cone having a sharp 



apex . The SFR in this galaxy is estimated to be in the range 2-5 M yr 1 (IFriedrich et al. 



2010|). Thus, the mechanical power from the stars is 1 — 4 x 10 42 erg s 1 , which is of the 



order of the kinetic power of the outflow. 



-43- 



A.4. NGC 3783 



An unresolved nuclear point source is detected in radio observations (Ulvestad & Wilson 
1984), ground-based observations in the optical (Winge et al. 1992), and HST [O ill] images 
presented in Schmitt et al. (2003). The total extent of the [O ill] emission is r ~ 0.6" (130 
pc). Our Br7 and [Si Vl] images also show an unresolved core but extended emission out 
to r = 140 pc along a PA^ 0°. In HST imag es, NGC 3783 presen ts the largest CLR, 



with more than 400 pc radius in the [Fe vil] line ( iMazzalay et al.ll2010l ). As pointed out by 
these authors, the compact [O ill] emission could be a result of an oversubtraction of the 
continuum image removing [O ill] emission from the outer regions of the halo. 

The Br7 morphology and velocity field are similar to that of [Si Vl], and rather different 
from the H2 1-0 S(l) (see Figure 113]) . The kinematics of the two regions are dominated 
by non-circular motions consistent with radial outflow. In this case the evidences are the 
following: (1) the radial velocity maps and the velocity channel maps show redshifted veloc- 
ities in the north part of the galaxy accelerating away from the nucleus out to ~ 0.35" and 
then appear to decelerate, (2) the maximum velocity of ~ 200 km s _1 observed at ~ 0.35" 
is too high to be explained by any reasonable gravitational potential and (3) the velocity 
dispersion of the ionized gas increases in the direction of the acceleration. A kinematic model 
incorporating rotation plus radial outflow provides a good fit to the data. The parameters of 
the model are presented in Table |5J The best fit PA disk = —15° and idisk = 38° are consistent 
with those of the H2 kinematics. The half-opening angle of the cone is 6*0^ = 34°, 6 m = 27° 
and its axis is inclined —30°. Therefore the northern cone is almost entirely redshifted and 
the southern cone is in its majority blueshifted, consistent with an intermediate viewing 
angle. Note that the blueshifted part of the bicone is in its majority hidden by the plane 
of the galaxy. The SW side of the galaxy is closer to us, and the southern cone therefore 
experiences more extinction. This is consistent with the Br7 and [Si Vl] images which show 
weaker, and in some place absent, emission SW of the location of the SMBH (Fig. [T5|) . The 
case of NGC 3783 provides strong evidence against the hypothesis that the radio jet is driving 
the bipolar outflow of i onized gas, since there is no collimated radio emission in the central 



parsecs of this galaxy (jOrienti fe Prietol 120091 ) so the clouds seem to accelerate even in the 



absence of any radio material. Furthermore, it is also the only Seyfert galaxy in our sample 
exhibiting broad lines and a line of sight far outside of the bicone (a = 26°). This gives 
support to a clumpy torus such that our view to the BLR is al most independent of the in- 



clinat ion and opening angles of the bicone. Within this context, IVollmer. Beckert fe Davies 



( 120081 ) found that NGC 3783 has the most transparent torus of the galaxies in our sample. 
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A.5. NGC 4151 



Pogge (1989) detected [0 ill] emission in ground-based images that is extended by 20" 
along PA= 50°. The [O ill] emission was later confirmed with HST to have a biconical 
morphology along PA= 55° (Evans et al. 1993). The total extent of the emission is 5.5" 
(350 pc). At radio wavelengths a complex structure of knots has been identified within the 
central few arcseconds along a PA of 77° (Mundell et al. 2003). 

Br7 and [Si vi] show a marginally resolved core as well as extended emission out to 80 pc 
along a PA=50°. Note that the complete extension of the NLR and CLR is not accesible in 
our observations due to the limited FOV of OSIRIS. The Br7 velocity field is similar to that 
of [Si vi], and rather different from the H 2 1-0 S(l). The kinematics of the two regions are 
dominated by non-circular motions consistent with radial outflow. In this case the evidences 
for outflow are the following: (1) the radial velocity maps and the velocity channel maps show 
redshifted velocities in the NE part of the galaxy and blueshifted velocities in the SW, which 
are accelerating away from the nucleus along a PA~ —130°, (2) the redshifted velocities 
in the NE follow the direction of the radio jet pointing away from us and the blueshifts 
in the SW are observed in the direction of the jet pointing towards us, (3) the maximum 
velocity of ~ 400 km s _1 observed at ~ 0.5" is too high to be explained by any reasonable 
gravitational potential and (4) the velocity dispersion of the ionized gas increases along a 
PA~ —130° consistent with the extended V-shape morphology observed in [O ill] images. A 
kinematic model incorporating rotation plus radial outflow provides a good fit to the data. 
The parameters of the model are presented in Table [51 The underlying rotation has a best 
fit PAdisk = 175° and idisk = 45°, which is consistent with the kinematics of the molecular 
gas in the nuclear region and roughly consistent with the large scale galactic rotation. The 
bicone has an external half-opening angle of 29°, an internal angle 9 m = 15° and is inclined 
^bicone = —45°. Therefore the NE cone is entirely redshifted and the SW cone is entirely 
blueshifted, consistent with the prediction of the Unified Model for a type 1.5 object. The 
misalignment of the NLR and the radio jet in NGC 4151 suggest that the radio jet is not 
the principal driving force on the outflowing NLR clouds in this galaxy. 



A.6. NGC 6814 



The [O ill] image of this galaxy (ISchmitt et al.l Il996l ) exhibits an unresolved core and 



some extended emission along PA~ —15°. The total extent of the emission is r = 0.6" (55 
pc). The Br7 and [Si vi] emission are consistent with this mophology, but in our OSIRIS 
data, the core is marginally resolved. This in fact is the most compact CLR in our sample. 
There is no emission related to the western radio extension found by Ulvestad & Wilson 
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(1984), but the overall distributions of Br7 and [Si vi] emission in the central ~ 100 pc seem 
to follow that of the radio emission (see Figure 121]) . 



The spatial resolution of the SINFONI data is estimated from the broad Br7 emission 
to be ~ 0.57" (FWHM), which is ~ 60 pc at the distance of NGC 6814. The azimuthal 
averages of the SINFONI observed distributions (Br7 and [Si vi]) have an HWHM of 30 
pc consistent with the PSF. The OSIRIS data has a higher spatial resolution of ~ 0.17" 
(FWHM), equivalent to 18 pc, and for this reason the properties measured from the OSIRIS 
data are taken to be more accurate, although, the measured and derived properties from 
both data sets are consistent, particularly the kinematics (see Figures l5l l6l [HI and [TO]) . Both, 
the Br7 and [Si vi] velocity fields, exhibit a perturbed rotational pattern. The best fit PA 
and i of the rotational component are consistent with those of the molecular gas kinematics. 
The velocity residuals map obtained after subtracting the rotational component to the [Si vi] 
(or Br7) velocity map reveal peculiar velocity gradients with blueshifts and redshifts on each 
side of the nucleus. These approaching and receding regions are also evident in the velocity 
channel maps (Figure [17]) . and indicate that material is accelerating away from the nucleus. 
Further evidence for outflows comes from the collimated radio emission which is spatially 
coincident with the coronal gas (see Figure |2"T1) . Assuming that the radio jet is approaching 
in the south part of the galaxy, the blueshifted velocities seen in [Si Vi] in this region would 
be outflowing. A kinematic model incorporating rotation plus radial outflow provides a good 
fit to the data. The parameters of the model are presented in Table [H Note that this model 
provides a good fit to both, the OSIRIS and the SINFONI datasets. The underlying rotation 
has a best fit PAd is k = 145° and idisk = 50°, which is consistent with the kinematics of the 
molecular gas in the nuclear region. The bicone has a half-opening angle 8 out of 47°, 6 m = 25° 
and is inclined bicone = —34°. Therefore the northern cone is almost entirely redshifted and 
the southern cone is almost entirely blueshifted, consistent with an intermediate viewing 
angle. Within the uncertainties, a direct view to the BLR is possible. 



A.7. NGC 7469 

A variety of observations at many wavelengths indicate that the active nucleus is sur- 
rounded by a r = 1.5" (480 pc) ring of starburst activity (Wilson et al. 1991; Mauder et 
al. 1994; Genzel et al. 1995; Soifer et al. 2003; Davies et al. 2004). The [O ill] emission is 
contained in a r = 2.5" region having an almost circular morphology. Within the SINFONI 
and OSIRIS FOV, which fall inside the ring, Br7 and highly ionized emission line gas, [Si Vi] 
and [Ca vm], are detected. In addition to the diffuse radio emission associated with the 
nuclear starburst ring, an E-W extended nuclear source is also detected (Thean et al. 2001), 
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which, at high spatial resolution, separates into five sources along an E-W direction within 
0.1" of the nucleus (Lonsdale et al. 2003). There is no indication in the Br7 flux distribution 
of a counterpart to this radio emission. However, the [Si Vl] emission seems to follow the 
direction of the nuclear radio sources. [Si vi] is extended along a PA=90° out to r = 90 
pc. The OSIRIS and SINFONI data are consistent with these measurements. Although the 
OSIRIS data have a slightly better spatial resolution than the SINFONI data of this galaxy, 
we decided to perform the kinematic analysis using the SINFONI datacube, because its FOV 
covers better the region where the CLR is extended (E-W). 

The Br7 emission detected in NGC 7469 is nearly symmetric, extends out to r = 125 
pc and peaks at ~ 0.1" (32 pc) west from the AGN location. [Si vi] presents a similar 
morphology, but it is slightly more extended in the E-W and peaks ~ 32 pc East from the 
AGN. These offsets may indicate that the Br7 emission is associated with star formation 
rather than the AGN. This is supported by the facts that the Br7 velocity dispersion is rather 
low in the central region and its velocity field is dominated by rotation. The velocity field of 
Br7 in the SINFONI and OSIRIS datasets is well fitted by a rotating disk with PA= 143° 
and % = 45°. These values are consistent with the PA and % of the molecular gas disk. The 
[Si Vi] velocity field, however, is rather different from the Br7 kinematics. It is dominated 
by a non-circular component exhibiting increments of velocity as a function of distance in 
the E-W direction. The magnitudes of the velocities reach a maximum at r ~ 0.16" and 
then appear to decrement. These are signatures of acceleration followed by deceleration and 
therefore indicate that the CLR is in an outflow. This phenomenon is observed in both, 
the OSIRIS and SINFONI datasets (see Figures 151 16} 151 and [TUl) . We performed a kinematic 
modeling of the CLR including rotation and radial outflow and found that the kinematics are 
best fitted with just an outflowing component. This can be clearly seen in the lower panels 
of Figure 1261 where we show the LOS velocities inside the bicone model as it would appear in 
the plane of the sky. The velocities outside this region do not show obvious regular patterns. 
The model parameters are presented in Table [5J The half-opening angle of the bicone is 
#out = 45°, the internal opening angle is &i n = 28 and its axis is inclined —33°. Therefore the 
eastern cone is almost entirely redshifted and the western cone is in its majority blueshifted, 
consistent with an intermediate viewing angle. As in the rest of intermediate type Seyferts 
of our sample (except NGC 3783), within the uncertainties, a direct view to the BLR is 
possible. 
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Fig. 1. — Example [Si vi] line profiles extracted from the central regions of three sample 
galaxies, chosen to represent the main types of profile shapes obseved over the fields. The 
data are shown by a solid black line, individual Gaussian fits by dashed blue lines and the 
continuum level fit in red. Below each spectrum is the residual plot. Fluxes and residuals 
are given in W m~ 2 /im" 1 . (a) An example line fit from NGC 6814 showing a symmetric 
profile. Similar profiles are observed over the fields of the majority of the galaxies (except 
NGC 1068); (b) a spectrum from NGC 3783 showing a slightly redshifted tail. These profiles 
are typical of NGC 3783 and NGC 6814; (c) example of a spectrum of NGC 2992 showing 
a slightly blueshifted tail. These profiles are typical of Circinus (see MS06), NGC 2992 and 
NGC 7469. 
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Fig. 2. — [Si vi] emission line profiles from different regions in the central 1.2" x 1.2" (80 x 80 
pc) of NGC 4151. The spectra have been extracted from apertures of 0.35" (23 pc) diameter. 
The locations of these regions are indicated by the arrows. The central panel shows the flux 
map of [Si Vi] emission. 
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Fig. 3. — [Si vi] emission line profiles from different regions in the central 4" x 4" (280 x 280 
pc) of NGC 1068. The spectra have been extracted from apertures of 0.5" (35 pc) diameter. 
The locations of these regions are indicated by the arrows. The central panel shows the flux 
map of [Si Vi] emission. 
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Fig. 4. — Two-dimensional maps of the Br7 flux distribution, velocity, and velocity dispersion 
(from left to right), for each galaxy observed with SINFONI. The galaxies are sorted by 
distance. The contours delineate the K— band continuum emission and the position of the 
AGN (the peak of non-stellar continuum emission at 2.2/im) is marked with a cross. Regions 
in white correspond to pixels where the line properties are uncertain and thus were masked 
out. These rejected pixels in the velocity and dispersion maps are those with a flux density 
lower than 5% of the peak of Br7 emission. 
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Fig. 5. — Continuation of Figure HI 
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Fig. 6. — Same as Figure HJ but for the galaxies observed with OSIRIS. 
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Fig. 7. — Azimuthally-averaged surface-brightness radial profiles of the Br7 flux distribu- 
tion for the galaxies observed with SINFONI and OSIRIS as indicated in the legend. The 
horizontal dashed line indicates the HWHMb^- Tipically, the standard deviation of the 
azimuthally averaged flux is 5%. 
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Fig. 8. — Two-dimensional maps of the [Si Vl] flux distribution, velocity, and velocity disper- 
sion (from left to right), for each galaxy observed with SINFONI. The galaxies are sorted by 
distance. The contours delineate the K— band continuum emission and the position of the 
AGN (the peak of non-stellar continuum emission at 2.2/im) is marked with a cross. Regions 
in white correspond to pixels where the line properties are uncertain and thus were masked 
out. These rejected pixels in the velocity and dispersion maps are those with a flux density 
lower than 5% of the peak of [Si Vl] emission. 
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Fig. 9. — Continuation of Figure [HI 
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Fig. 10. — Same as Figure El but for the galaxies observed with OSIRIS. 
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Fig. 11. — Flux maps of [Si vi], [Al ix] and [Ca vm] emission (second, third and fourth 
columns) in Circinus, NGC 1068 and NGC 2992. The first column represents the continuum 
emission at 2.2/im. These are the only three galaxies in our sample exhibiting these three 
near-IR coronal lines. Note the similarities in their morphologies and that [Al ix] presents 
always the most compact emission. See text for details. 
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Fig. 12. — Azimuthally-averaged surface-brightness radial profiles of the [Si vi] flux distri- 
bution for the galaxies exhibiting [Si vi] emission. The horizontal dashed line indicates the 
HWHM[si V il- Tipically, the standard deviation of the azimuthally averaged flux is 5%. 
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Fig. 13. — Velocity fields of H 2 (left), Br7 (middle) and [Si Vl] (right) emission in the central 
regions of NGC 3783 (FOV~ 160 pc), NGC 4151 (FOV~ 160 pc), NGC 6814 (FOV~ 250 
pc) and NGC 7469 (FOV~ 250 pc). The [Si Vl] velocity fields show increments of velocity 
with distance in a different direction to the H 2 rotation indicative of an outflow of > 100 km 
s _1 . See text for details. 
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Fig. 14. — Velocity tomography of the [Si vi] emission-line in the Seyfert 2 galaxy NGC 
1068. The channel maps were obtained by integrating the flux within velocity bins of 200 
km s _1 along the [Si vi] emission-line profile. The number in the upper part of each panel 
correspond to the central velocity of the bin in km s _1 relative to systemic. The bottom 
right panel shows the integrated intensity image by summing the velocity slices. 
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Fig. 15. — Velocity tomography of the [Si Vl] emission-line in the Seyfert 1.5 galaxy NGC 
3783. The channel maps were obtained by integrating the flux within velocity bins of 80 
km s _1 along the [Si vi] emission-line profile. The number in the upper part of each panel 
correspond to the central velocity of the bin in km s _1 relative to systemic. The bottom 
right panel shows the integrated intensity image by summing the velocity slices. Note that 
the acceleration in this galaxy is mostly observed in the north part. See text for details. 
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Fig. 16. — Velocity tomography of the [Si vi] emission-line in the Seyfert 1.5 galaxy NGC 
4151. The channel maps were obtained by integrating the flux within velocity bins of 120 
km s _1 along the [Si vi] emission-line profile. The number in the upper part of each panel 
correspond to the central velocity of the bin in km s _1 relative to systemic. The bottom 
right panel shows the integrated intensity image by summing the velocity slices. 
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Fig. 17. — Velocity tomography of the [Si vi] emission-line in the Seyfert 1.5 galaxy NGC 
6814. The channel maps were obtained by integrating the flux within velocity bins of 80 
km s _1 along the [Si vi] emission-line profile. The number in the upper part of each panel 
correspond to the central velocity of the bin in km s _1 relative to systemic. The bottom 
right panel shows the integrated intensity image by summing the velocity slices. 
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Fig. 18. — Velocity tomography of the [Si vi] emission-line in the Seyfert 1.5 galaxy NGC 
7469. The channel maps were obtained by integrating the flux within velocity bins of 80 
km s _1 along the [Si vi] emission-line profile. The number in the upper part of each panel 
correspond to the central velocity of the bin in km s _1 relative to systemic. The bottom 
right panel shows the integrated intensity image by summing the velocity slices. 
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Fig. 19. — Velocity tomography of the [Si vi] emission-line in the Seyfert 1.9 galaxy NGC 
2992. The channel maps were obtained by integrating the flux within velocity bins of 80 
km s _1 along the [Si vi] emission-line profile. The number in the upper part of each panel 
correspond to the central velocity of the bin in km s _1 relative to systemic. The bottom 
right panel shows the integrated intensity image by summing the velocity slices. 
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Fig. 20. — Velocity tomography of the Br7 emission-line in the Seyfert 2 galaxy Circinus. 
The channel maps were obtained by integrating the flux within velocity bins of 80 km s" 1 
along the Br7 emission-line profile. The number in the upper part of each panel correspond 
to the central velocity of the bin in km s _1 relative to systemic. The bottom right panel 
shows the integrated intensity image by summing the velocity slices. 
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Fig. 21. — Contour plots of radio continuum emission overlaid on the [Si vi] flux maps. The 
overlays correspond, from left to right, to NGC 1068 (radio image from Gallimore et al. 
1996), NGC 4151 (radio image from Mundell et al. 2003) and NGC 6814 (radio image from 
Ulvestad & Wilson 1984). The legend "blue" indicates the direction of the radio jet coming 
towards us as well as the regions where the blueshifts in the velocities of the coronal gas are 
observed. 
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Fig. 22. — Example of a bicone model incorporating radial acceleration and deceleration for 
two test galaxies having different inclinations and PAbi CO nc = 0°. Left: Three-dimensional 
structure of the model with the arrow indicating the LOS. Top Right: Front and back 
projections of a galaxy having ibicone = 90°. Bottom right: Front and back projections of 
a galaxy having bicone = 0°. The velocities are assigned a colour based on the amount of 
redshift (warm colours) or blueshift (cold colours). Green represents approximately zero 
velocity. 
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Fig. 23. — Maps showing the data and the best-fit kinematic model used to reproduce the 
[Si Vl] kinematics in NGC 3783. The maps are centered at the position of the AGN (black 
cross), and the scales are in km s _1 . Top Left: [Si Vl] velocity field used for the modeling, Top 
Middle: Best-fit kinematic model incorporating rotation and biconical outflow, Top Right: 
Residuals (data- model) to the fit, Bottom Left: Rotational component of the combined 
kinematic model shown in the top middle panel, Bottom Right: Outflowing component of 
the combined kinematic model shown in the top middle panel. 
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Fig. 24.— Same as Figure [23] but for NGC 4151. 
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Fig. 25.— Same as Figure [23] but for NGC 6814. 
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Fig. 26. — Maps showing the data and the best-fit kinematic model used to reproduce the 
[Si Vl] kinematics in NGC 7469. The maps are centered at the position of the AGN (black 
cross), and the scales are in km s _1 . While the top row shows the LOS velocities in the 
entire FOV, the bottom row shows only the LOS velocities inside the bicone model. Left 
column: [Si Vl] velocity field used for the modeling, Middle column: Best-fit kinematic model 
corresponding only to biconical outflow, Right column: Residuals (data-model) to the fit. 
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Fig. 27. — Maps showing the data and the best-fit with pure rotation for a rotation- 
dominated object (NGC 6814, bottom), and the best fit with pure outflow for an outflow- 
dominated case (NGC 3783, top). The maps are centered at the position of the AGN (black 
cross), and the scales are in km s -1 . Left column: [Si vi] velocity fields used for the mod- 
eling, Middle column: Best-fit incorporating only one kinematical component (top-outflow, 
bottom- rotation), Right column: Residuals (data- model) to the fit. In the two cases, the 
quality of the fit decreases when compared to the best fits incorporating rotation plus bicon- 
ical outflow. 
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Fig. 28. — Maps showing the data and the best-fit kinematic model used to reproduce the 
[Si Vl] kinematics in NGC 2992. The maps are centered at the position of the AGN (black 
cross), and the scales are in km s _1 . Top Left: [Si Vl] velocity field used for the modeling, Top 
Middle: Best-fit kinematic model incorporating rotation and biconical outflow, Top Right: 
Residuals (data- model) to the fit, Bottom Left: Rotational component of the combined 
kinematic model shown in the top middle panel, Bottom Right: Outflowing component of 
the combined kinematic model shown in the top middle panel. 
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Fig. 29. — Graph showing the anticorrelation found between the half opening angle of the 
outflow (0 O ut) and its maximum velocity (v max ). The error bars correspond to the uncertain- 
ties in the estimation of the best-fit parameters (Tables [5] and [6]). 
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Fig. 30. — Graph showing the anticorrelation found between the half opening angle of the 
outflow (0 O ut) an d the molecular gas mass (M H2 ) in the central r < 30 pc. The values of M H2 
were obtained from H09 (except for NGC 2992 which was taken from Friedrich et al. 2010). 
The vertical error bars represent the typical 35% error of the gas mass estimates (H09). The 
horizontal error bars correspond to the uncertainties in the estimation of # out (Tables [5] and 
ED- 
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Fig. 31. — Graph showing the correlation found between the maximum velocity of the outflow 
(ii max ) and the molecular gas mass (Mh 2 ) in the central r < 30 pc. The values of Mh 2 were 
obtained from H09 (except for NGC 2992 which was taken from Friedrich et al. 2010). The 
vertical error bars represent the typical 35% error of the gas mass estimates (H09). The 
horizontal error bars correspond to the uncertainties in the estimation of v max (Tables [5] and 

ED. 
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Fig. 32. — Geometric models of the NLR/CLR bicones and the inner galactic disks in NGC 
1068 (top left), NGC 2992 (top right), NGC 3783 (middle left), NGC 4151 (middle right), 
NGC 6814 (bottom left) and NGC 7469 (bottom right) based on the parameters listed in 
Tables |5] and [6j The color blue indicates that the majority of the radial velocities in the cone 
are blueshifted. Accordingly, a red cone indicates that the majority of the radial velocities 
are redshifted. See Section 15.41 for details. 
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Table 1: Summary of AGN data included in this study 



00 
05 



a Classification is taken from the catalogue of IVeron-Cettv fc Veronl ( 20061 ). and Distance (D) from the 

NASA/IPAC Extragalactic Database. 
6 Spatial resolution (FWHM) estimated from the data itself, using the broad Br7 emission and the non-stellar 

continuum, except for Circinus (see Section 2.2). The numbers in parentheses indicate the FWHM in pc. 
c Column indicating if the line was Detected (D), Marginally Detected (MD) or Not Detected (ND) in the 

spectra. 
d The FOV in these observations has a cross shape. 
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Table 2: Flux, size and velocity dispersion of the NLR 



a Nuclear Br7 flux integrated in a circular aperture of 2xFWHM (two times the resolution) and centered at 

the nucleus. The numbers in parentheses indicate the diameter of the aperture in pc. 

Total Br7 flux integrated in a rectangular aperture covering the observed emission (inside 2xRbi- 7 ). 
c Photometric semimajor axis of Br7 emission measured using as reference the contours corresponding to 5% 

of the peak of Br7 flux. Lower limits are indicated by a colon. 
^Photometric semimajor axis of [O in] emission 
e References for R[ p iii] and PA[q ,;;]: (1 ) Marconi et al. (1994), (2) ISchmitt et al.l (|2003l ). (3) Garcia-Lorenzo 

et al. (2001), (4) lschmitt et all |l996J ). (5) This work. 
•^The PA in this galaxy changes from ~ 10° in the central 50 pc to ~ 110° at larger scales. 
3 Measured in the 640 x 640 pc SINFONI datacube (see Fig. 27 of D07.) 
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''Size and PA of the compact bright core. The radius of the starformation ring is ~ 1300 pc. 



Note. — Uncertainities in the flux measurements are approximately 5% 
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Table 3: Flux, size and velocity dispersion of the CLR 



"Nuclear [Si vi] flux integrated in a circular aperture of 2xFWHM (two times the resolution) and centered 
at the nucleus. The numbers in parentheses indicate the diameter of the aperture in pc. 

fc Total [Si vi] flux integrated in a rectangular aperture covering the observed emission (inside 2xR[gj v ;j). 

c Photometric scmimajor axis of [Si Vi] emission measured using as reference the contours corresponding to 
5% of the peak of [Si vi] flux. Lower limits are indicated by a colon. 

^Measured in the 640 x 640 pc SINFONI datacube (see Fig. 27 of D07.) 



Note. — Uncertainities in the flux measurements are approximately 5% 
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Table 4: Best-fit kinematic parameters of a rotating disk 



a Mean value of the absolute values in the velocity residuals map. /i > 20 km s l indicates a bad fit. 
^Values obtained from H09, except for NGC 2992 ( Friedrich et al.luOlOT ). These values are fou n d to be 
consistent with those of the large scale galactic disk. This was verified by H09 and iFriedrich et al.l (|2010l ). 
c Values affected by artifacts in the SINFONI detector (see Figure HJ). 

^Results obtained keeping PA and i constant during the fit, alt hough the data indic ate a change in the PA 
from i-v 0° in the central r < 0.7" to ~ 40° in the outer region ( Friedrich et al.ll2010l) . 



Note. — The typical estimated error of the kinematic parameters is ±12 c 
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Table 5: Results of kinematic modeling of the CLR in five Seyfert galaxies. 



a In this galaxy the fitting procedure incorporates only radial outflow and is described in a separate publication 

(Miiller- Sanchez et al. in prep.). 
& No deceleration observed. 

c Full module of the outflow velocity vector, i.e. relative to the nucleus not to the observer. 
d Radial distance from the AGN to the location where f max is observed. 
e Rotational velocity measured at r — 50 pc. 
•^Angle between the outer edge of the bicone and the LOS. 
^Relative tilt of the bicone axis to the disk. In NGC 1068 and NGC 7469 this angle is estimated using PAh 2 

and in 2 (see Table @|. 
''Minimum x 2 values for models incorporating only rotation. 
'Minimum \ 2 values for models incorporating only biconical outflow. 
■'Minimum \ 2 values for models incorporating rotation and biconical outflow. 
fc Number of data points in the velocity fields, equivalent to the number of degrees of freedom. 



Note. — The quoted parameters correspond to the model presenting the smallest (x 2 ) ■ and the uncer- 
tainties refer to 3c confidence level. The derived values are also the best-fit parameters for the NLR, except 
in NGC 7469. 
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Table 6: Kinematic model of the CLR in NGC 2992 

a No deceleration observed. 

b Full module of the velocity vector, i.e. relative to the nucleus not to the observer. 

c Radial distance from the AGN to the location where f max is observed. 

^Rotational velocity measured at r = 50 pc. 

e Angle between the outer edge of the bicone and the LOS. 

^Relative tilt of the bicone axis to the disk. 

9 Minimum \ 2 value for models incorporating only rotation. 

h Minimum \ 2 value for models incorporating only biconical outflow. 

^Minimum \ 2 value for models incorporating rotation and biconical outflow. 

^Number of data points in the velocity field, equivalent to the number of degrees of freedom. 



Note. — The quoted parameters correspond to the model presenting the smallest (% 2 ) . and the uncer- 
tainties refer to ia confidence level. The derived values are also the best-fit parameters for the NLR. 
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Table 7: Summary of results of the kinematic analysis 



"Column indicating if the corresponding kinematical component is Dominant (D), Present (P) or Not Present 
(NP) in the velocity field. "U" stands for uncertain. The asterisk indicates that in these galaxies the size 
of the NLR, as traced by [O ill] emission, is larger than our FOVs. However, the model parameters of 
NGC 1068, NGC 2992 and NGC 4151 are consistent with those derived in previous kinematic studies of the 
NLR in these galaxies with larger FOVs. All CLRs are fully covered by the SINFONI and OSIRIS FOVs, 
except in the Circinus galaxy. See text for details. 

^Column indicating the table in the paper in which the best fit kinematic parameter values are presented. 
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Table 8: Mass outflow rates (M out ) and their corresponding kinetic powers (E ont ) compared 
to mass accretion rates (M acc ) and accretion powers (L bo \) 



a ateral surface area of the biconical outflow obtained directly from the kinematic models, using r t and # ut- 
References for L hol : (1) Prieto et al. (2010), (2) IWoo fc Uml (J2002I ). 



